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Abstract 


This thesis is devoted to properties of quark matter at high density and (compar¬ 
atively) low temperature. In nature matter under these conditions can possibly 
be found in the interior of (some) neutron stars. Under "normal" conditions 
quarks are confined in the hadrons. Under the extreme conditions in the core 
of a neutron star, however, the strong interaction between the quarks becomes 
weaker because of asymptotic freedom. Then the hadrons will break up, and for 
the remaining interaction between the quarks a (semi-)perturbative treatment 
should be sufficient, at least at asymptotic densities. 

It turns out that cold dense quark matter does not behave like a Fermi liq¬ 
uid as a consequence of long-range chromomagnetic interactions between the 
quarks. While the specific heat of a Fermi liquid is a linear function of the tem¬ 
perature (T) at low temperature, the specific heat of normal (non-superconduc¬ 
ting) quark matter is proportional to T log T at low temperature. In this thesis 
the specific heat and the quark self energy in normal quark matter are discussed 
in detail. A discrepancy between earlier papers on the specific heat is resolved 
by showing that in one of these papers gluonic contributions to the energy den¬ 
sity were overlooked. Moreover higher order corrections to the known leading 
order terms are computed in this thesis in a systematic manner. 

The quark-quark interaction mediated by one-gluon exchange is attractive in 
the color antitriplet channel. Therefore quarks form Cooper pairs at sufficiently 
small temperatures, which leads to the phenomenon of color superconductivity. 
During the last years the issue of gauge (in)dependence of the color supercon¬ 
ductivity gap has been discussed in the literature. Here a formal proof is given 
that the fermionic quasiparticle dispersion relations in a color superconductor 
are gauge independent. In a color superconductor the gluon field acquires a 
non-vanishing expectation value in general, which acts as an effective chemical 
potential for the color charge. This expectation value is computed for two dif¬ 
ferent color superconducting phases (2SC and CFL). It is shown that at leading 
order this expectation value is determined from a tadpole diagram with a quark 
loop. Tadpole diagrams with gluon or Nambu-Goldstone loops turn out to be 
negligible at our order of accuracy. 
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"Die Astronomen miissen aus dem Strahlen den Stern erkennen - (und tap- 
pen doch im Dunkeln herum, denn sie konnen durch ihr Verfahren, man nennt 
es Spectralanalyse, doch nur die erdverwandten Stoffe auffinden - das dem Stern 
Ureigene bleibt ewig unerforschlich) ..." 


- Gustav Mahler in a letter to Alma Schindler 




Chapter 1 

Introduction 


We know four fundamental interactions in nature, namely gravitational, electro¬ 
magnetic, weak and strong interactions. Since this thesis is devoted to proper¬ 
ties of quark matter, we shall mainly be concerned with the strong interactions. 
Today the accepted theory of strong interactions is Quantum Chromodynamics 
(QCD). It is a gauge theory with gauge group SU (3) c , where the index c stands 
for color. The quarks belong to the fundamental representation of this group, 
and the gauge bosons (gluons) belong to the adjoint representation. Denoting 
the quark field with ip and the gluon field with A, one writes the Lagrangian of 
QCD as [1, 2] 

1 Nf 

Cqcd = - -F« V F^ + (i#ij - mfSij) Aj, ( 1 . 1 ) 

/=i 

with 0 = 'fDp, F“ v = d, x Al ~ d v A% - gf abc A b ^A c v: and D, M] = d^j + igA^T*. 

In Eq. (1.1) we have written explicitly the sum over quark flavors. Apart from 
the invariance under SU( 3) c gauge transformations, the Lagrangian (1.1) is also 
invariant under global unitary transformations in flavor space as long as the 
differences of the quark masses can be neglected. In the chiral limit (mf —► 0) the 
Lagrangian is even invariant under separate flavor space transformations of left 
and right handed quarks, with the corresponding symmetry groups SU(Nj)l 
and SU(Nf)n. Moreover the Lagrangian (1.1) is invariant under global U( 1) 
phase transformations of the quark fields, which is related to quark number 
conservation 1 . 

The quantum theory based on the Lagrangian (1.1) is renormalizable [4,5, 6]. 
Neglecting the quark masses for the moment, one finds for the one-loop beta 
function 2 

_/%)_= (U - !"/) ■ (1-2) 

1 At the classical level the Lagrangian is also invariant under axial phase transformations of the 
quark fields, which constitute the group 17(1) a - In the quantum theory, however, this symmetry 
is broken by an anomaly, see e.g. [3]. 

2 By now the beta function is known to four-loop order [ 7 ], 


1 



Nature chooses Nf = 6. (For many physical situations the number of "active" 
flavors is smaller than six. E.g. for the astrophysical systems, on which we 
eventually focus, we need to take into account only the three lightest flavors.) 
Therefore the beta function is negative, which shows that QCD possesses the 
all-important property of asymptotic freedom [8, 9]. 

Experimental confirmations of QCD come from deep inelastic scattering ex¬ 
periments, and other high energy processes (see [10] and references therein). 
These experiments give a current world average of [10] 

a s {M z ) = 0.1187 ± 0.002 (1.3) 

with a s = |^, and M z = 91.1876 ± 0.0021 GeV is the mass of the Z boson. 

At small energies quarks and gluons are confined into hadrons. In the con¬ 
fined phase chiral symmetry is spontaneously broken by a non-vanishing quark 
condensate (qq), with 7r, K and rj playing the role of (pseudo-)Nambu-Goldstone 
bosons. The low-energy dynamics of the hadrons can be described with chiral 
perturbation theory [11,12]. 

From asymptotic freedom one expects at sufficiently high temperatures a 
phase transition to a deconfined phase [13]. The resulting "soup" of quarks 
and gluons is called quark gluon plasma (QGP). This idea is corroborated by 
lattice simulations, which predict a phase transition at T c ~ 170 MeV for Nf = 2 
[14,15], and at T c ~ 150 MeV for Nf = 3 [15]. It is assumed that the QGP existed 
in the early universe. Experimentalists try to reproduce the QGP in heavy ion 
collisions at SPS (CERN), at RHIC (Brookhaven), and from 2007 onwards at LHC 
(CERN). 

The equation of state (the pressure as a function of the temperature) of the 
high temperature QGP has been computed perturbatively up to order g 6 log g 
(see [16] and references therein 3 ). The expansion in the coupling constant turns 
out to be only poorly convergent for non-asymptotic temperatures. It is however 
possible to improve the situation by performing systematic resummations of 
the perturbative series, such as resummations within dimensional reduction [16, 
19], HTL-screened perturbation theory [20, 21], or approximately selfconsistent 
resummations based on the 2PI effective action, [22, 23, 24] which show good 
agreement with the lattice data for T > 2 T c . For reviews on the high temperature 
QGP see e. g. [25, 26,27, 28, 29,30]. 

It is of great importance to understand the properties of the quark gluon 
plasma also at finite chemical potential 4 . This region of the QCD phase diagram 

3 The coefficient of the <? 6 -term cannot be computed within perturbation theory, since Feynman 
diagrams of arbitrarily high order contribute to this coefficient [17,18], and there is no known way 
of how to resum them. 

4 The chemical potential can be viewed as a measure for the density at a given temperature. 
E.g. for an ideal relativistic gas of free fermions the particle number density is proportional to /i 3 
for T —> 0. Thus a high chemical potential is equivalent to high particle number density in this 
case. 
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is relevant for the interior of compact stars [31, 32, 33] 5 . From asymptotic free¬ 
dom one expects a phase transition to a deconfined phase to occur not only at 
high temperature, but also at high density. Unfortunately one does not know 
an efficient algorithm for lattice simulations at finite //. The reason is that the 
fermion determinant in the path integral is complex at finite /j, which precludes 
Monte Carlo importance sampling. So far this "sign problem" can only be cir¬ 
cumvented for n < T (see [28] and references therein). But for the region // 3> T, 
which is relevant for compact stars, no method of implementing lattice simu¬ 
lations is known to date. Therefore simplified models, such as Nambu-Jona- 
Lasinio models [34, 35, 36], and (semi-)perturbative methods are the only viable 
approaches at the moment. In compact stars the chemical potential might be as 
large as 500 MeV. At this energy scale a s (n) is of the order one, therefore the ap¬ 
plicability of (semi-)perturbative methods is rather questionable. Nevertheless 
one can try to extrapolate results which have been obtained for small values of 
the coupling constant to larger values of the coupling constant, hoping that the 
qualitative features remain valid. 

Compared to the chemical potential, the temperature in the interior of a com¬ 
pact star will be rather small (of the order of tens of keV [32]). In this region of 
the QCD phase diagram quark matter is expected to be in a color supercon¬ 
ducting phase, which is characterized by a non-vanishing diquark condensate. 
Such diquark condensates can arise from the fact that one gluon exchange is 
attractive in the color antitriplet channel, leading to the formation of Cooper 
pairs at sufficiently small temperatures. Color superconductivity has been dis¬ 
cussed already in the late 1970's [37, 38], but at that time it was believed that 
the color superconductivity gap would be of the order 1 MeV, and thus almost 
negligible. Only some years ago [39, 40] it was realized that the gap may be as 
large as ~ 10 — 100 MeV, with corresponding critical temperatures of the order 
~ 5 — 50 MeV. This discovery stimulated extensive research in this field during 
the last few years, see e.g. [41, 42, 43, 44, 45, 46, 47, 48] for reviews. 

Due to the interplay of finite quark masses and the constraints from color 
and electric neutrality, there may be various color superconducting phases in 
the QCD phase diagram, which are distinguished mainly by the particular form 
of the diquark condensate. These matters will be reviewed in more detail in 
chapter 5 of this thesis. 

In nature color superconducting phases could in principle be discovered 
from analyses of neutron star data. There exist indeed some proposals for sig¬ 
natures that could indicate the presence of color superconductivity. For in¬ 
stance, the cooling behavior of a neutron star depends on the specific heat and 
the neutrino emissivity, which are both sensitive to the phase structure (see 
e.g. [49, 50, 51, 52, 53]). Other possible signatures include r-mode instabili- 

5 As far as the macroscopic structure of compact stars is concerned, the curvature of spacetime 
is in general not negligible. However, the change in the metric is tiny for the typical length scales 
of particle physics [31]. Therefore it is sufficient to consider the field equations of matter in flat 
spacetime. 


3 



ties [54, 55], and pulsar glitches due to crystalline structures in inhomogeneous 
color superconducting phases [56]. 

This thesis is organized as follows. In chapter 2 we discuss some formal 
properties of gauge theories, which will be used in chapter 5. In chapters 3 and 
4 we compute the quark self energy and the specific heat of normal degenerate 
quark matter. These results may be relevant for the cooling properties of (proto-) 
neutron stars [57], Chapter 5 is devoted to color superconductivity. We give a 
general proof that the fermionic quasiparticle dispersion relations in a color su¬ 
perconductor are gauge independent. Furthermore we compute gluon tadpole 
diagrams, which are related to color neutrality. Chapter 6 finally contains our 
conclusions. 

Throughout this thesis we shall use the following conventions. We use nat¬ 
ural units, h = c = ks = 1. The Minkowski metric is g ,iu = diag(l, —1, —1, —1). 
Four momenta are denoted as K 1 ' = (fc 0 ,k). The absolute value of the three 
momentum is denoted as k := |k|. A unit three vector is denoted as k := k /k. 
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Chapter 2 


Symmetries in quantum field 
theory 

2.1 General gauge dependence identities 

In this section we will recapitulate the derivation of a general gauge dependence 
identity, following [58, 59, 60]. 

Let us consider an arbitrary gauge theory which is defined by an action func¬ 
tional Si nv [<p] that is invariant under some gauge transformations, 

v = ^]<r- (2.i) 

Here we use the DeWitt notation [61], where an index i comprises all discrete 
and continuous field labels, and a Greek index (a) comprises group and space- 
time indices. E.g. for QCD, as defined by the Lagrangian (1.1), one has explicitly 
[with a <-> ( x, a)] 

D l a -»• ( 5 ab dl - gf abc A c ^(y))5(x -y), i <-► A^y), 

-igT a ip{y) 5{x -y), V ’(y), (2-2) 

ig${y)T a S(x -y), ip(y). 

In general we assume that the gauge generators form an off-shell closed al¬ 
gebra, which means that the commutator of two gauge transformations is again 
a gauge transformation, 

(2-3) 

where the "structure constants" cG can in principle be field dependent. 

In order to quantize the theory we have to fix the gauge freedom, e.g. with a 
quadratic gauge fixing term, 

SinvYp\ S g f[<p\ = S inv [ip\ + ^F a [ip\F a [p]. (2.4) 
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The quantum theory is defined via the path integral representation of the gener¬ 
ating functional of connected Green's functions 1 , 

exp [iW[J]) = j Vpdet (F?[p\D l p [p\) exp (* (S gf [<p\ + Jup 1 )) . (2.5) 

With the help of a Legendre transformation 

W[J\=T[p\ + J i p\ (2.6) 

one obtains the effective action, 

exp(iT[<p]) = j Vp det (F?[p\D l p[p]) exp (* (S gf [p] - T^[p\(p - (p) 1 )) . (2.7) 

If a change 5F a [<p] in the gauge fixing condition is accompanied by the fol¬ 
lowing (non-local) gauge transformation, 

V = D l a [p\8Z a [p\, = g a pWF%], (2.8) 

the gauge fixed action S g f will remain invariant. Here G°g is the ghost propaga¬ 
tor in a background field tp, which is defined via 

M</>] = -8«. (2.9) 

Let us examine in which way the measure and the Faddeev Popov determinant 
change under the gauge transformation (2.8) and 5F a . From the measure we get 
a Jacobian, which can be evaluated for small 8p l using the formula 

det(l + X) ~ 1 + TrX, (2.10) 

which is valid if \X t j\ <C 1. The variation of the Faddeev Popov determinant can 
be calculated using Eq. (C.l). In this way we find 

<5 (Vp det (F?£>y [p]) = Vp det (FfD},) Tr^^ + «5Tr log (F^D},) 

= Vpdet (F“^) [d^SF? + gOpF? {D^D^ - D^D^J) G\5F 5 ] . (2.11) 

The first term between the brackets in the second line vanishes if I)' ry i = 0. Using 
Eqs. (2.3) and (2.9) one finds that the remaining part vanishes provided that 
Cp a = 0. These two conditions are fulfilled for QCD 2 , and also for the Higgs 
models which we will discuss in this chapter. 

: We assume for the moment that p 1 are bosonic fields in order to keep the notation simple. It 
is easy to check, however, that the final identity Eq. (2.12) is also valid for the fermions in QCD, 
provided one chooses suitable conventions for the fermionic derivatives. 

2 In order to prove this one uses the fact that the structure constants are totally antisymmetric 
and that the Gell-Mann matrices are traceless. 
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The only term in Eq. (2.7) which is not invariant is the source term 
Therefore we find the following gauge dependence identity, 

sm = =: ( 2 . 12 ) 

where we use the notation 

( K W\)W\ 

= e _irM J T>(pK[p] det {F?[(p\D t p[<p]) exp (i (, S gf [ip] - T *[<p](<p - <p) 1 )) (2.13) 
for any K [p ]. 

As an application let us consider the effective potential V e ff for a translation- 
ally invariant system, 

r[p]|^ = const. = ~Veff{(p) J d A X. (2.14) 

In general V e ff will be a gauge dependent quantity. However, the identity (2.12) 

ensures that 

SVeff{0) + 5X'[y] dVe ^) = o. (2.15) 

This is the Nielsen identity for the effective potential [59, 62], which states that 
a change in the gauge fixing function (the first term on the left hand side) can 
be compensated by a change in <p (the second term on the left hand side). This 
identity implies that the value of the effective potential at its minimum (where 
dV e jf / dtp 1 — 0) is gauge independent. 

A further consequence of the identity (2.12) is the gauge independence of the 
position of the propagator singularities [60, 62, 63, 64], which we will discuss in 
more detail in the sections 2.3 and 2.4 and in chapter 5. 

2.2 Implications of global symmetries 

The aim of this section is to examine the consequences of global symmetries for 
the one-point and two-point functions in quantum field theory. 

Let us consider an arbitrary quantum field theory with field content ( V(x), 
T'fx'). We assume that the theory is invariant with respect to a group G of con¬ 
tinuous global symmetry transformations, which act on the fields as 

&i(x) = (2.16) 

where the g t] form some representation of G. The effective action up to second 
order can be written as 

r[<J>] ~ J d 4 xTi(x)®i(x) + J d 4 xTi(x)$i(x) 

+ f AAr^x.^iW^!/) + 0(4> 3 ). (2.17) 


7 



We assume that the symmetry is not (spontaneously) broken. Then we have 
T[$] = T [$'], which implies 


r i{x) = gjiTj(x), (2.18) 

r i (x) = r j (x)gl j , (2.19) 

Fij(x, y) = gkiT k i(x, y)g] v (2.20) 

for all g. First let us consider Eqs. (2.18) and (2.19). With the additional as¬ 
sumption that the representation formed by the set of all g's is non-trivial and 
irreducible, these equations imply that 


Ti(s) = 0, fi(x) = 0. (2.21) 

which means that the tadpole diagrams vanish (in other words, the expectation 
values of the field operators are zero). 

If we assume that the representation is unitary, we can write Eq. (2.20) as 

T ji(x, V) = gJ^Vf k {x, y)g ki ■ (2.22) 

If we assume furthermore that the representation is irreducible, we can now 
invoke Schur's lemma, which gives the result 

T,J oc Sij, (2.23) 

i.e. the (inverse) propagator is proportional to the unit matrix with respect to the 

group indices. 

Next we would like to examine in which way parity (P) and time reversal 
symmetry (T) constrain the propagator. For simplicity we work in Euclidean 
space and assume that our fields <p l transform under PT as 3 

<?(x) -»• cWtfi-x), (2.24) 

with cW = ±1. [This comprises for instance vector fields and (pseudo-)scalar 
fields.] Assuming translational invariance, the bilinear part of the effective ac¬ 
tion can be written as 

m\ = \f ^<t>\-K)T l3 (KW{K), (2.25) 

where T t] is the full inverse propagator. Invariance under PT implies 

m\ = j ^VWTijiKWi-K). (2.26) 

3 No summation over the index in parentheses. 
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Performing the substitution K —► —it' in the last equation and comparing it with 
Eq. (2.25) we find r ij(K) = c (>> e <J) T, : j(—K). The mere definition of the inverse 
propagator implies Tij(K) = Tj^—K), and therefore we have 

Tij{K) = Sc^TjiiK). (2.27) 

Using Dyson's equation 

A'-'P,,, = -4 (2.28) 

we arrive at the following identity for the full propagator, 

S ij (K) = c (i) c (i) S ji (I \). (2.29) 

For c^cCO = l the propagator is therefore symmetric in i and j. 

2.3 Abelian Higgs model 

This and the next section serve as an illustration of the application of the gauge 
dependence identities for systems with spontaneous symmetry breaking at fi¬ 
nite temperature. We will discuss the gauge independence of the locations of 
propagator singularities, both for the Abelian and a non-Abelian Higgs model. 
This is actually more than a simple exercise, since the methods that we develop 
here will be useful also in the more complicated context of color superconduc¬ 
tivity, which will be the subject of the last chapter of this thesis. 

Let us also mention that the static limit of the gauged Higgs model yields the 
Landau-Ginzburg Lagrangian, which can be used for instance for the descrip¬ 
tion of superconducting systems in the vicinity of the transition temperature 
[38]. In this thesis, however, we shall not go into the details of the Landau- 
Ginzburg description of superconducting systems. 

There exists a vast literature on various aspects of Higgs models, let us just 
mention Refs. [65, 66, 67, 68, 69, 70, 71, 72], in which the finite temperature 
behavior of Higgs models is discussed. 

2.3.1 Definition of the model 

The Abelian Higgs is defined by the Lagrangian 

Qnv = -\F^F^ + - U($) 

with F /im = d^A u - d v A^, D fl = td^ + ieA p/ <f> = ^(<£i + 2 ), and 

D($) = 

The Lagrangian (2.30) is invariant under the gauge transformation 

(5d> = ia&, 5A^ = — d^a. 


(2.30) 

(2.31) 

(2.32) 
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We assume m 2 < 0 so that the gauge symmetry is spontaneously broken. The 
expectation value of $ is taken to be real so that is the Higgs boson and y >2 
is the would-be Nambu-Goldstone boson. C inv is invariant with respect to a 7l 2 
symmetry [59]: 

Afj, —> -Ap, p 2 ->• -p '2 (2.33) 

Following [59] we choose a gauge fixing which does not break the 7L 2 symmetry: 

C gf = ■ A +K(p 2 ) 2 (2.34) 

with some constant k and / /t = <9 /t , f fl = ( 0 . <),) or /), = u /t , where u 7i = 
(1,0, 0,0). In momentum space we write / /t = /3(K)K^ + 7 (Ji)h M with /V' = 
( 5 ^ _ K»K v /K 2 )u„. If ffi = Ufi we will assume k = 0 because otherwise time 
reversal symmetry would be violated. The corresponding ghost Lagrangian is 
given by 

£ g host = c^--^f-d + 2ampi\c. (2.35) 

The total Lagrangian 

-d = C-inv T •£<;/ T t~-ghost. (2.36) 

is still invariant under (2.33), which implies that Green functions (derivatives of 
the effective action, T |a=^ 2 =o,^i=^) with an odd total number of external A- 
and <p 2 -legs will vanish. Using Dyson's equation (2.28) is is easy to see that this 
is also true for general Green functions. 


2.3.2 Free propagators 

The free propagators can be obtained easily by evaluating the path integral for 
the free theory. One finds that the free propagators are simply given by inverting 
the "coefficients" of the kinetic (bilinear) terms in the Lagrangian. If (pjfx) is the 
solution of the free field equation. 


SSba. 

5(fi(x) 


= Mx), 


one finds for the free propagator 


Gij(x,y ) = 


Hj(y) 

5Ji(x) 


(2.37) 


(2.38) 


In this way one obtains the following tree level propagators for the Abelian 
Higgs model (with the notation "5"=tp 2 ): 


i_ 

(0) K 2 - e 2 (p 2 


( A _ 

v /3iv 2 — enip 2 0) 


^(o) 


C> JV 
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_ 7 2 fe 2 + {aK 2 + k 2 )(K 2 - e 2 ^ Q) ) , 

(f3K 2 - eKip (0) ) 2 )' 

ik 

K 2 (I< 2 - e 2 tf 0) )m 2 ~ e^ ( o)) 2 
x [ (-e 7 2 /c 2 <^ (0 ) - K 2 ((3k + ae(p [0) ){K 2 - e 2 <pj 


(2.39) 


A (0) 


+— (2 (3 2 K 2 k + a/3eK 2 (f( 0 ) + ae 2 K<p (0 )) , 

e 2 7 2 A: 2 ^ 0) + J\ 2 (K 2 - e 2 <p 2 {0) )((3 2 K 2 + ae 2 cp 2 {0) ) 
K 2 (K 2 - e 2 cpf 0) ){(5K 2 - emp (_ 0 )) 2 


K 2 ~ m 2 - |(^ 2 0) ’ 


where 


A^ = g^ - 


K ll K u n^n v 


B I,M = 


1 K 11 K v 

C' 1V = —(h fJ ‘K u + W), V IJV = 

k K z 

£^=—, = 


and <p(o) is the position of the minimum of the tree level potential, i.e. 

/—6m 2 


(2.40) 


(2.41) 


(2.42) 


(2.43) 

(2.44) 

(2.45) 


(2.46) 


We notice that all the structure functions are gauge dependent at tree level apart 
from A a(0) , A b(0) and Ajg*’ 1 . 

The (full) propagators have the following symmetries: the definition of prop¬ 
agators implies A y = A 7 *, and therefore in the non-condensed notation 

A • U/ (K) = A^(-K), A f*(K) = A 5/i (—AT). (2.47) 

Invariance with respect to parity and time reversal implies (see Sec. 2.2) 

A ^(K) = A^(-K), A ^(K) = -A^i-K). (2.48) 

Eqs. (2.47) and (2.48) yield 

A '"'(AT) = A^(K), A ^(K) = -A^ l {K). (2.49) 

If one used the gauge condition (2.34) with / /t = u /t and n / 0, (2.48) and (2.49) 
would not be valid, and an explicit calculation shows that a term which is anti¬ 
symmetric in fiv would indeed appear already in the tree level propagator A^. 
We may parametrize the full propagators as follows: 


A**" = Aa^ + AbB^ + AcC^ + AdP' 11 ', 

A m5 = Ai^ + As^, 


(2.50) 

(2.51) 
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We perform an analogous decomposition for the two-point functions T^ and 
F /( 5 which are obtained from the effective action. These are one-particle irre¬ 
ducible apart from possible tadpole insertions. We also define self energies in 
the usual manner as IIy = — A r. 1 . 

With the help of Dyson's equation (2.28) one can express the structure func¬ 
tions of the full propagators in terms of the self energies in the following way: 


A a 

A B 
Ac 
A D 
Ai 
A 2 
A55 


with 


Ta’ 

_ rf - r 55 rch 2 _ 

-r?r B + 2T 1 t 2 t c + r 2 r D + r 55 (r 2 , + r B r D )h 2 ’ 

Tir 2 + r 55 rch 2 

r 2 - r 55 r D h 2 s ’ 
r| + r 55 r B ra 2 . 
rf - r 55 r D h 2 
(TiTb - r 2 r c )h 2 
r 2 - r 55 r D h 2 
(TiTg - r 2 r D )h 2 
r 2 -r 55 r D h 2 s ’ 

(r 2 c + r B r D )n 2 A 

r 2 - r 55 r D h 2 


(2.52) 

(2.53) 

(2.54) 

(2.55) 

(2.56) 

(2.57) 

(2.58) 


Ta = 

-(K 2 -e 2 tf 0) ) + U A , 

(2.59) 

r B = 

I^!_ (i y2_ e 2-2 o)) + nB) 

(2.60) 

r G = 

^+n c , 

a 

(2.61) 

r D = 

P 2 X 2 , 2 _ 2 
a + e P(0) + n ^’ 

(2.62) 

id = 

/ K0\ , „ 

-* 1 ep(0) + — 1 fc + ni, 

(2.63) 

r 2 = 

-^h 2 + n 2 , 

a 

(2.64) 

h55 = 

2 ft 2 

A' 2 + — + n 55 . 

(2.65) 


2.3.3 Gauge independence of propagator singularities 

Taking the second derivative of the gauge dependence identity (2.12) one ob¬ 
tains 

ST,ij = T Mj 8 X k + r M 5 X k + T kj i 8 X k + T k SX ^. (2.66) 
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We evaluate (2.66) at A = p 2 = 0, p\ = p. Then the last term vanishes. For the 
inverse Higgs propagator we obtain in configuration space 


s 2 r 


5<pi(x)5(pi(y) 


= dz 


<s 3 r 


8ipi(z)6(pi(x)5<pi(y) 


_5X^ ( z ) 


+ 


f dz _Al 


ir^-\ 5XV1 W + (*<-*/)• W 

(x)\<pd<pi(y) <p 


Using translation invariance the first term on the right hand side can be rewrit¬ 
ten as 


dz 


d' 3 r 


5ipi(z)5pi(x)5pi(y) 


0 ^( 0 ) 


d 


5 2 T 


v dip 6ipi(x)8ipi(y) 


SXw ( 0 ) 


In momentum space we obtain therefore 


(TT tp li p 1 + Sp — 2T ipiipi 6X%\ 


( 2 . 68 ) 


(2.69) 


or 


+ S ^ A X, = 


(2.70) 


with Sip := —dX p] (z = 0). In a similar way one obtains for the gauge field 
propagator 

SA^ + Sp-^A pu = -Aw8X" p - A up 8X p - A p5 SX" 5 - A v5 6X£, (2.71) 


which implies 


d 


5A-/ + 51P—A-/ = A 


(2.72) 


(5A7 1 + Sip^A-t = 2A 


a 


-1 


dp 




h u h p rir 2 + r 55 r c n v K p \ 

~ O I -^.O „ „ _ O , ) u u\ , 


n* 


r? - r 55 r D 


n 


+ 


Fir c - TaT^ 

rf - r 55 r D h 2 


n v 6X% 


(2.73) 


Eqs. (2.70), (2.73) and (2.72) imply that the locations of the poles of the Higgs 
propagator and of the transverse and longitudinal components of the gauge 
field propagator are gauge independent [62], provided that the singularities of 
the SX's do not coincide with those of A vitpi , A a or A b, respectively. In the case 
of A^ one also has to take into account the various F's and h's on the right hand 
side of (2.73) [64, 60]. As in [64, 60] one may argue that 8X is 1PI up to a full 
ghost propagator and tadpole insertions. The singularities of the ghost propa¬ 
gator will be different from the physical dispersion laws since there is no HTL 
ghost self energy. The tadpoles are 1PI up to a full Higgs propagator evaluated 
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at zero momentum, and at this point the Higgs propagator is non-singular. In 
general, the 1PI parts of SX are expected to have no singularities, apart from pos¬ 
sible mass-shell singularities which can be avoided by introducing an infrared 
cut-off to be lifted only at the very end of the calculation [73, 74], 

On the right hand side of (2.73) additional singularities arise from the factor 
1 / K 2 contained in the h's and from 1 /(T 2 — TssT nn 2 )- These kinematical singu¬ 
larities have to be excluded from the gauge independence proof. The expression 
T j — r 55 r i/n 2 is obviously gauge dependent already at tree level, and in general 
according to 


S[rj - r 55 r D h 2 } + s<p-^[ r? - r 55 r D h 2 ] = 


= 2 


r 2 — r 55 r D h 


2\K ,l Kp , r _ r „\ KI ' n P \SX p p 


K 2 - (n r 5 5rc + rir 2 ) : 


+ ((riTc - ^D^2)n p — 2T 1 T£) 


frK, 


6X p 5 - (Tj + n 2 T 55 T D )5X% . (2.74) 


2.4 A non-Abelian Higgs model 

Let us consider a non-Abelian Higgs model with the gauge field and the scalar 
field both belonging to the adjoint representation of SU(2). The Lagrangian is 
given by 

Cinv = --F^aK" + (D^) a (D»<p) a - V(<p), (2.75) 

with Fa U = d p A v a - fXA'a + ge ahc A p A u c , (D p ^) a = d p (p a + geabcA^c, and 

V(<p) = + ^(<Pa<Pa) 2 , (2.76) 

with m 2 < 0. We assume that the expectation value of <p a is proportional to 5 a 3 . 
Using the discrete symmetries of the Lagrangian [59], it is easy to see that the 
gauge field propagator is diagonal in color space. It should be noted, however, 
that the gauge field propagator in the Higgs phase is not proportional to the unit 
matrix, which can be expected from the discussion in Sec. 2.2. 

It turns out that the Higgs propagator X ip3lf3 fulfills the same gauge depen¬ 
dence identity as in the Abelian case [Eq. (2.70)]. For the gauge field propagator 
one finds that and both fulfill Eqs. (2.71)-(2.73), and that A 33 obeys the 
same gauge dependence identity as in QED [60]. The fact that the gauge depen¬ 
dence identities are identical with those of Abelian models is of course only due 
to the particular simplicity of the SU ( 2 ) model. 

The SU( 2 ) Higgs model is an interesting toy model, since we will see in 
chapter 5 that also in the case of color superconductivity the gluon propagator 
is in general not proportional to the unit matrix in color space. 
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Chapter 3 

Quark self energy in 
ultradegenerate QCD 

3.1 General remarks 

In this chapter we will compute the quark self energy in normal ultradegenerate 
QCD, both at zero and (small) finite temperature. As discussed in the introduc¬ 
tion, and in more detail in chapter 5, quark matter at high density and suffi¬ 
ciently small temperature is actually in a color superconducting phase. Never¬ 
theless, it is interesting to study also the properties of normal, i.e. non-supercon- 
ducting quark matter. First, a good understanding of normal quark matter is cer¬ 
tainly helpful before tackling the more complicated case of color superconduc¬ 
tivity. Moreover, it is conceivable that quark matter in young (proto-)neutron 
stars is already dense enough for deconfinement, but still not cold enough for 
color superconductivity. If this is the case, the results of our computations will 
be relevant in particular for the cooling behavior of proto-neutron stars. 

The quark self energy will be used in the next chapter to compute the specific 
heat of normal cold dense quark matter. We will evaluate the quark self energy 
only on the light cone, since this is the quantity which enters the formula for the 
specific heat at leading order [see Eq. (4.15) below]. Furthermore it is sufficient 
for the computation of the specific heat to consider the fermionic quasiparticles 
in the vicinity of the Fermi surface. Thus the quark momentum is hard, while 
the momentum of the gluon in the one-loop self energy (see Fig. 3.2) may be 
arbitrarily soft. Therefore [75] we have to use a resummed gluon propagator in 
order to avoid an IR divergence in the quark self energy. For this reason we will 
review the computation of the HDF resummed gluon propagator in the next 
section. 
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Figure 3.1: One-loop gluon self energy. 


3.2 Gluon self energy 


The gluon self energy at the one-loop level is given by the diagrams of Fig. 3.1. 
The second one of the four diagrams depends on the chemical potential and 
the temperature, while the three remaining ones give only finite temperature 
corrections to the gluon self energy Below we will only need the gluon self 
energy at zero temperature, therefore we evaluate here only the diagram with 
a quark loop in Fig. 3.1. For QCD the diagram contains an additional factor of 
Tr [T a T b ] = rjd ah compared to the photon self energy in QED. In the following 
we shall drop the Kronecker delta, and define 


9eS = < 


9 2 N f 

2 


9 2 N f , 


QCD, 

QED. 


(3.1) 


Then our final results will be applicable both for QCD and QED. 

In the imaginary time formalism [76] we get for the second diagram of Fig. 
3.1 

v [■ p3 u 

IW*w, q) = - 9 2 eS T £ j ^3 Tr [^S 0 (K - Q) lv S 0 {K)] , (3.2) 


where Q^ 
by 


(ico, q) and K 1 ' = (iuj n . k), and the free quark propagator is given 


So(iu) n i 



dk 0 ftpojK) 

2ir ko — iuj n — p 


(3.3) 


with 

Po(K) = 27rsgn (k 0 ) 6 (k% - k 2 ) = ^ [5(k 0 - k) - S(k 0 + k)}. (3.4) 

Using the methods outlined in appendix A it is straightforward to evaluate the 
Matsubara sum in Eq. (3.2), 




1 1 

k 0 - iio n — ppQ — iio n + iu> - p 


nfjkp - p) - nf(po - p) 
ko-Po- iu 


(3.5) 


where nj is the Fermi-Dirac distribution defined in Eq. (A.13). After performing 
the Dirac trace in Eq. (3.2) we get [77], 





^-po(ko,k)p 0 {p 0 ,p) 
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(3.6) 


x ^ - n ’ {K - '*> ^ ■ p ), 

— 4*0 — iu} 

where p = k — q. Now we can perform the analytic continuation itv — > qo + +) 
to get the retarded self energy. Using 


-J— = V--8(x) 
x + irj x 


(3.7) 


we can separate the real and imaginary parts of the self energy. It is convenient 
to define 

n G := cT n M „, Uh := uWILfu,, (3.8) 

with a 1 ' = (1,0,0, 0). These two structure functions are related to the transverse 
and longitudinal self energies via 


n L =-^-yr-njT, n T = 1 (-n L + n G ). 

q- 2 


(3.9) 


After subtraction of the vacuum contribution one finds for the real parts [78] 


ReU G (qo,q) = / dkn(k) 

7 r z 


x ( 4 k + 


> o 

?0 “ <? 2 | (2 k + q 0 + q)(2k - q 0 + q) 


2q 


■log 


(2k + q 0 ~ q)(2k - q 0 ~ q) 


ReU H (q 0 ,q) = [ dkn(k) 

n 2 Jo 

(2k + q 0 + q)(2k + go - q) 
4 q 

(2k - q 0 - q)(2k -qo + q) 


2k ( 1 — — log 

qo + q 

L V q 

1 

o 


+ - 


4 q 


log 


log 


2k + q 0 + q 


2k + qo~q 
2k- qo~q 


2k -q 0 + q 


where we use the abbreviation 


1 


n(k) = - (n f (k - n) + n f (k + /z)). 


(3.10) 


(3.11) 


(3.12) 


For the imaginary parts one finds [77] 


imn G, H (qo,q) = °’ g ) “ n G,H(~Po,p)), 


with 


(3.13) 


7r G (9o ,q) = —l dkn(k)(ql - q 2 )F(k,q 0 ,q), (3.14) 

7T p Jo 

^n(qo,q) = [ dkn(k)(2k + q-qo)(2k-q-q 0 )F(k,q 0 ,q), (3. 

2vrp J o 


15) 
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where 


F(k,q 0 ,q) := sgn(k - qo)Q(\k - q\ < \k - q 0 \ < \k + q\), (3.16) 

where we use the notation ©(true expression) = 1 and ©(false expression) = 0. 

In general the remaining integrations over k can be performed analytically 
only for the imaginary parts [77]. In the zero temperature limit, however, a 
closed result can be obtained also for the real parts. In this limit one finds for for 
II l and n T for go > 0 [77] 


Relit (go, q)\ T=0 = [ 32 ^ + ^(90, q ) 


ReII T (go,g)| T=0 = 


9eS 


16p 2 g (q 2 + 2ql) 


967 r 2 q 3 

+ {q 2 ~ ql) (6g 2 M(g 0 , q) ~ N(q 0 ,q )) 

g 2 s {q 2 - ql) 


Imn L (go,g)| T =o = ~ 


Imn T (g 0 ,g )| r=0 = 


487rg 3 

W(q 0 ,q)U(q 0 ,q) -W(-q 0 ,q)U(-q 0 ,q) 

5eff ( q 2 - ql) 


96irq 3 


{W (g 0 , g) + 6g 2 ) U (g 0 , g) - (W (~q 0 , g) + 6g 2 ) U (-g 0 , g) 

where we have used the abbreviations [77] 


(3.17) 


(3.18) 


(3.19) 


(3.20) 


U(qo,q) = (2p- |g 0 + g|)0(2p - |g 0 + g|), (3.21) 

R(go,g) = (2p-g-go)(2p + 2g-g 0 ), (3.22) 

IR (go, g) = @(go + g)R(go,g) + ©(-go-g)R(-go,-g), (3.23) 

R(qo,q) = (go + g)log|go + g| - (2p + g 0 + g)log|2p + g 0 + g|, (3.24) 

S(qo,q) = (—12p 2 g 0 + (2g — g 0 )(g + g 0 ) 2 ) log |g 0 + g| 

+(2 p -2g + go)(2p + g + go) 2 log |2p + g + g 0 |, (3.25) 

M{qo,q) = R(qo,q) - R(qo,-q) +R(-qo,q) - R(-qo,-q), (3.26) 

N(qo,q) = S(q 0 ,q) - S(q Q ,-q) + S(-q 0 ,q) - S(-q 0 ,-q). (3.27) 


If g, go <C p the leading contribution to the self energy is quadratic in //. This 
contribution comes from hard momenta in the quark loop (k ~ //), therefore 
the corresponding self energy ist called Hard Dense Loop (HDL) self energy 
[75, 79, 80, 81]. In the following we shall denote the HDL approximation with a 
tilde ('). The explicit expressions are considerably simpler than the general ones 
given above [76], 


Refl L (go, g) 


2m 


2 r 


«0-fUog 
2q 


qo + q 


go 


(3.28) 
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Ren T (qo,q) 


Imfli (go, g) 
lmU T (q 0 ,q) 


2 <?o ( qo , q 2 - ql , 

™ ~ I — + log 


m 


q \q 2 q 2 

; qo (q 2 - ^ 


qo + q 


I qo-q 

2 „2 




7T0(g- - g 0 ), 


- «§) 


2g d 


(3.29) 

(3.30) 

(3.31) 


where m 2 = g 2 s p 2 /(2 tt 2 ). The leading finite temperature corrections have the 
same functional dependence on go an d q, one only has to include a finite tem¬ 
perature correction into m 2 . For QCD one finds 


2 r g 2 g 2 
m 2 = + 

J An 2 




(3.32) 


where the contribution proportional to N c comes from the gluon and ghost loops 
in Fig. 3.1. The part of the gluon self energy which is quadratic in the tempera¬ 
ture is called Hard Thermal Loop (HTL) self energy. 

Using Dyson's equation II pu = D~J — D^j iu one can obtain the one-loop 
resummed gluon propagator. In the HDL approximation one finds in covariant 
gauge [76], using the tensors defined in Eqs. (2.43) and (2.44), 


D/aJqo- q) — DtA^u D/, // 




9 9 

q 2 - qo 


and in Coulomb gauge (see e.g. [82]) 


_ _ q2, _ g2 

D /w (qo, q) = -DtA^, u + D l A— — g^ 0 g u0 + 


where 


D l = 


Dx = 


q- 


1 


q 2 -ql a n L 

1 


q 2 - ql + n T 

The HDL spectral density for the transverse part is given by 

Pr{qo,q) = 2lmD T (qo + ir],q). 


(3.33) 

(3.34) 

(3.35) 

(3.36) 

(3.37) 


In covariant gauge the longitudinal propagator contains an unphysical factor of 
l/(g 2 — ql) [60, 64], therefore one defines the longitudinal spectral density as 

2 2 

PLiqo, q) = 2 9 ° q ImDfJqo + ip, q). (3.38) 

q 2 

Both px and pi, contain a pole and a cut contribution. The explicit expressions 
can be found in [76], and are also given in appendix B for convenience. In terms 
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Figure 3.2: One-loop quark self energy with resummed gluon propagator. 


of the spectral densities, the structure functions of the propagator can be written 
as [76] 


Dr(qo, q) 

q 


f°° dq' 0 p T (q' 0 ,q ) 

7-oo 27T q' 0 -qo 
f°° dqopUqjq) _ J_ 
7-oo 27T q'o -qo q 2 ' 


(3.39) 

(3.40) 


3.3 Fermion self energy on the light cone 

The fermion self energy is defined through 

S-\P) = Sf\P) + X(P), (3.41) 

where Sq(P) = is the free fermion propagator, and F’ 1 ' = (p°, p) with 

p° = iuj n + p,ui n = (2n + 1)ttT in the imaginary time formalism, and F’ 1 ' = ( E , p) 
after analytic continuation to Minkowski space. 

With the energy projection operators 

A± = l(l± 7o7 y) (3.42) 

we decompose £(P) in the quasiparticle and antiquasiparticle self energy, 

£(P) = 7oA+S + (P) - 7o A-S_(P), (3.43) 

and 

yoS- 1 = S+ X A+ + SZ l A- (3.44) 

so that Sf 1 = — [p° T (|p| + £±)]. 

The one-loop fermion self energy is given by 

£(P) = - 9 2 C f Tj2 f -0^rS o {P-Q)YD, v {Q), (3.45) 

where D /u/ is the gauge boson propagator. 

Following [83] we introduce an intermediate scale q* , such that rri <C q* <C //, 
and we divide the ^-integration into a soft part (q < q*) and a hard part (q > q*), 

S + = EJ ) + Sf ) . (3.46) 
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For the hard part we can use the free gluon propagator, whereas for the soft part 
we shall use the HDL resummed gluon propagator. 

The hard part of the fermion self energy on the light cone in a general Cou¬ 
lomb gauge is given by the gauge independent expression 

E ( ± ) (E) = — J dqq 2 j dt J dr 0 [5(r 0 - r) - 6(r 0 + r)] 

X / 2 (sgn(r 0 ) — p-q rq) J dq 0 ^[6(q 0 - q) - S(q 0 + <?)] 

1 + n b {q 0 ) - n f (r 0 - p) - 9(g 0 ) + 0(-r o ) 
r 0 + qo ~ E - ir) 

+ (sgn(r 0 ) + r-p) [™/(^o ~ d) ~ @(-?’o)] |, (3.47) 


where r = p — q and E = p. The distribution functions rtf, and rif are given 
by Eqs. (A.5) and (A.13). In Eq. (3.47) we have subtracted the vacuum parts 
of the distribution functions, since we know anyway that the vacuum part of 
the fermion self energy on the light cone vanishes because of gauge and Lorentz 
invariance. This subtraction is in fact necessary, because otherwise we would 
get a spurious vacuum contribution coming from the fact the we do not use an 
0(4) invariant cutoff for the energy-momentum integration. The finite-T and 
finite-// parts, however, are not affected by these subtleties in the UV region. 

After performing the angular integration we find 


Th) 


{E ) = f ]t dqq^2(q-\p-q\+pcoth(^) 

(l -)_ e (|p-g|-/d/ T ) (i + e (b-d+M)/ T ) 


+ 7d°g | ^ e ( p+ q_ M )/T 1 j ^ e (p + g +M )/Tj 

The integral is finite in the limit q* —> 0, with the result 

xV h) = ^f(i + o(T 2 /p 2 )), 


(3.48) 


(3.49) 


with = g 2 Cfp 2 /(A-K 2 ). Hereg* enters as a correction proportional to (g*/p) 2 , 
so that we can send g* to zero. Conversely we expect that in the soft contribution 
we should be able to send q* to infinity without encountering divergences, as 
will indeed be the case, but only after all soft contributions are added together. 

For the soft part one finds on the light cone in a general Coulomb gauge the 
gauge independent expression [82] 


S { 1\e) 


J dqq 2 J dt J dr 0 [<5(r 0 - r) - S(r 0 + r)] 
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x< 2 (±sgn(?’ 0 ) - p-qr-q) 


+ (±sgn(?’ 0 ) + r-p) 


r ^ PT ( qo , q ) i+nb{qo) ~ n f {r °-^ 

-oo 27T 

dqo , , 

pL{qo,q ) 


r 0 + ?o =F \E\ - ill 
1 + n b (q 0 ) - n/(r 0 - p) 
ro + qoE \E\ - irj 


1 1 


-n f (r 0 -p) 


(3.50) 


where E = ± p, and pt and pl are the spectral densities of transverse and longi¬ 
tudinal gauge bosons, as given in Eqs. (3.37) and (3.38). We may use q « \E\,r 
because of q < q* and \E\ ~ p. Depending on the sign of E, we can drop the 
term S(r o + r ) or the term 5(r o — r) in Eq. (3.50), since its contribution is sup¬ 
pressed with ~ q/E compared to the remaining contribution. Then we find for 
the soft contribution to the real part of £ + 


Rp 5]( s ) — — ^ 

Ke + “ 8vr 2 


dqq 2 


dt 


/-1 


— [(1 - t 2 )p T {qo,q) + PL(qo,q)\ 

IT 

(3.51) 


X -p 1 +"<■(«>-"/(fi-M-gt) _ 1 (1 _ 2 {E _ _ qt)) 

qo - qt <r 


This quantity vanishes for E = p by symmetric integration. After performing 
the ^"integration we therefore have 


KeL + - 47r 2 


dqq 2 


i-i 


dt ( rif(E — p — qt) — rif(—qt)) 


x(l -t 2 ) \ReD T {qt,q)-ReD L {qt,q) 


(3.52) 


where Dt and Dl are given in Eqs. (3.35) and (3.36), respectively. For InT 
(which receives no hard contribution) we find in an analogous way 


Im£ + = — 


s 2 c f 

87T 2 


dqq 2 J ^ dt ((1 - t 2 )p T {qt , q) + pL(qt, q)) 


x [1 + n b (qt) - rif(E - p - qt)] 


(3.53) 


(s) 

The antiquasiparticle self energy £_ is obtained by inserting negative values 

(s) 

of E in the above expressions for T . and including an overall factor (—1). With 
p > 0 we can then replace rif(E — p — qt) by 1. 


3.4 Expansion for small \E — fi\ and small T 

In this section we will perform an expansion of £ + in the region 

T ~ \E - p\ <C gp <C p. (3.54) 

This region is relevant for the computation of the low temperature specific heat, 
see chapter 4. We will use the expansion parameter a := T/m, and we define 
A := (E — p)/T. From (3.54) we have a<l and A ~ 0(1). 
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3.4.1 The first few terms in the series 

In the part with the transverse gluon propagator we substitute 

q = ma 1/3 z, t = a 2/3 v/z. (3.55) 


After expanding the integrand with respect to a we find for the transverse con¬ 
tribution 


ReS 


(s) _ g 2 C f ma 


. Q 
am 


+(T) 


dv 


IRA 


dz 


e A — 1 


7 T- 


Ja 2 / 3 \v\ 

am. 


2v 2 z(v 2 it 2 — 4 z b ) o/o 
+ —^ a 2/3 


\_v 2 ir 2 + Az & (v 2 tt 2 + 4z 6 ) 2 


(1 + e v )(l + e x ~ v ) 

16r; 4 z 3 (3t> 2 7r 2 — 4 z 6 ) 
(v 2 tt 2 + 4 z 6 ) 3 


a 4 / 3 + ... 


(3.56) 


The ^-integrations are straightforward. In the u-integrals we may send the inte¬ 
gration limits to Too. Using the formulae 1 


dv 


e x — 1 
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(1 + e v )(l + e x ~ v ) 1 
= T(a + 1) Li a _(_i(—e A ) — Li a+ i(—e A ) 
e A — 1 
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■ log |u| 
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(3.57) 


(3.58) 


we find, neglecting terms which are suppressed at least with ( m/q *) A , 


= — g 2 Cfm. 

[ a 


x - 


+ 


12vr 2 
2 i/ 3 a 5/s 


A log ( 2 ^ q J- 1 +7^A- ^ ^Li a+ i(-e A ) — Li a+ i(—e A )^ 
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9^3vr 7 /3 
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+ 


8(24 — 7r 2 )a 3 logo 
27^ 


A(A 2 + 7T 2 ) + 0{a 3 ) 


(3.59) 


1 r(a) is the Euler Gamma function, defined as T(a) := The Gamma function 

satisfies T(n) = (n — 1)! for n € N and T'(l) = — 7 e, where 7 e is Euler's constant, 7 e = 
lirrim^oo (EfcL 1 1 ~ logm) ~ 0.5772156 ... Furthermore, Li a (s) is the polylogarithm, defined as 

Li a (s) = Er=U fc fc-“. 
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In the longitudinal part we substitute q = nix and t = au/x. In a similar way 
as for the transverse part we find 


= —g 2 Cfm 


aX 
87 t 2 


log 
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(q 


* \2 


(7T 2 — 4)a 3 log a 
967T 2 


A(A 2 + 7T 2 ) + 0(a 3 ) 


(3.60) 

Turning now to ImS + we notice that it vanishes at E = // only in the case of 
T = 0. For finite temperature, however small, there is an IR divergent contribu¬ 
tion in the transverse sector [76], 


T v, , 9 2 C f T , m 

ImE +(T) | E= „=-^-log — 


(3.61) 


where the infrared cutoff may be provided at finite temperature by the nonper- 
turbative magnetic screening mass of QCD. In QED, where no magnetostatic 
screening is possible, a resummation of these singularities leads to nonexponen¬ 
tial damping behavior [84], 

After subtraction of the energy independent part we have 


ImX 1 — ImS_ 
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i£i 

87r 2 


dqq 2 / dt (nf(E —/j, — qt) — nf(-qt)) 
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x [(1 - t 2 )PT(qt, q) + PL(qt, q)] ■ 


(3.62) 


Following the steps which led to Eq. (3.56), we find for the transverse contribu¬ 
tion 
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(3.63) 


Using the formula 2 


dv- 


e A — 1 


r— |v|“sgn(a) 


(1 + e v )(l + e x ~ v ) 

= — r(a + 1) [”Lio,_|_i(—e A ) + Li a _|_i(—e A ) + 2(1 — 2 Q ) £(a + 1) 


Vet > 0 

(3.64) 


we find in a similar way as above 

ImS+ ( D - ImS +( T )lfi=/i = ^C/mj - Y^ lo g cosh Q) 

2 C(o.) is the Riemann zeta function, which is defined as Q(oi) := k~ a = Li a (l). 
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- W r ® + Li V3(- A ) + 2 (l - f ©) 

- 20 ^W r (I) (Li7/s(-e- A ) + Li 7/ 3<-^) +2(1- 2- 4/3 ) C (D) 

+0(o 3 )l. (3.65) 


For the longitudinal part we obtain 

r 2 \2 

ImS +(L) - ImS^ L) | B=A( = -g 2 C f m + 0(a 3 ) . (3.66) 

Putting the pieces together, and using the abbreviation e = E — p, we obtain 
for the real part 



where 


/t(A) 

1 d 

= exp 1 - +-—^Li Q:+ i(-e” A ) - Li Q+ i(-e A )^j 

1 , (3.68) 

O !=0 

/ 2 (A) 

= r ( 1 ) ( Li 5 / 3 (-e _A ) - Li 5 / 3 (-e A )) ' , 

(3.69) 

/s(A) 

= -T (t) (Lh/3( -e A ) - Li 7 / 3 (—e A )^ , 

(3.70) 

MX) 

= |A(A 2 + 7 t 2 )| 1/3 . 

(3.71) 


The determination of the function 03 (A) requires resummation of IR enhanced 
contributions, which will be discussed in the next subsection. We note further¬ 
more that the dependence on q* indeed drops out in the sum of the transverse 
and longitudinal parts. 

The functions /*(A) show a simple asymptotic behavior. In the zero temper¬ 
ature limit (| A| —> oo ) 3 we have /,(A) —> |A|. If the temperature is much higher 

3 This limit can be taken by replacing the distribution functions in Eq. (3.52) with their zero 
temperature counterparts, or by using in Eqs. (3.68)-(3.71) the asymptotic expansion of the poly¬ 
logarithm given e.g. in [85]. 
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than |E — p| (i.e. A —> 0) we have /i(A) —> cq := |exp(l - 'je) = 2.397357... 
and / 2 , 3 , 4 (A) —> 0. For |A| S> Co or |A| <C Co we may approximate /i(A) with 
max(co, |A|), which is qualitatively the result quoted in [ 86 ]. It should be noted, 
however, that the calculation of Ref. [ 86 ] only took into account transverse gauge 
bosons, and therefore the scale under the logarithm and its parametric depen¬ 
dence on the coupling was not correctly rendered. 

For the imaginary part we find 


ImS_|_ - ImS + | E = g 2 Cfm. 
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(3.72) 


where 



51(A) = 2 log cosh Qj , 

(3.73) 

32(A) = 

(|) ( Li 5/s(“ e A ) + Li5/3( _eA ) + 2 (l - 2 C (§)) 

1 3 /5 

, (3.74) 

33(A) = 

A|, 

(3.75) 

34(A) = 

(t) ( Li 7/3( _e A ) + Li7/3(-e A ) + 2 ^1 - 2 4// ' 3 ^ ( (|) 

V1 3/7 

)J -(3-76) 


In the zero temperature limit we have 5 , (A) —> |A|. If the temperature is much 
higher than | E — // we have g, (A) —> 0. Again the determination of the func¬ 
tion (^(A) requires resummation of IR enhanced contributions, which will be 
discussed in the next subsection. 

Explicitly, our T = 0 result reads 
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(3.77) 


where we have included the zero temperature limits of C 3 and ci.-j, which will 
be derived in the next subsection. We observe that the self energy at T = 0 is 
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at leading order proportional to e log \e\, which is nonanalytic at e = 0. This be¬ 
havior is a consequence of quasistatic chromomagnetic fields, which are only dy¬ 
namically screened. On the other hand chromoelectric fields, which are screened 
by the Debye mass, give an analytic term (oc e) in the self energy at leading order, 
and nonanalytic contributions (oc e 3 log |e|) only at higher orders. 

Apart from the first logarithmic term, the leading imaginary parts contri¬ 
buted by the transverse and longitudinal gauge bosons were known previously 
[82, 87, 88 ]. As our results show, the damping rate obtained by adding these 
two leading terms [82, 87] is actually incomplete beyond the leading term, be¬ 
cause the subleading transverse term of order |e | 5 / 3 is larger than the leading 
contribution from Eg,. 


3.4.2 Evaluation of the a 3 -coefficient 


We write the function 03 (A) in Eq. (3.67) as the sum of the transverse and the 
longitudinal contribution, C 3 = C 37 + c^l- The terms through order a 3 in the 
integrand of Eq. (3.56) give the following contribution to C 3 t(A), 
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^A(A 2 + 7T 2 ) 

^4 + Ji (" 8 ( n - 61 °g (§)) + 27f2 (! - lQ g (§))) 


8(24 - 7 r 2 ) 
97T 6 


: (3 — 2yij)A(A 2 + 7r 2 ) 



(3.78) 


Evaluating explicitly the next few terms in the expansion of the integrand, one 
finds additional contributions of order a 3 . They arise from the fact that the 
^-integrations would be IR divergent, being screened only by z mm = a 2 ^\v\. 
Therefore also terms which are formally of higher order in the integrand con¬ 
tribute to the order o 3 in the self energy. In the following we will perform a 
systematic summation of all these terms. 

With the substitution (3.55) the transverse gluon self energy in the HDL ap¬ 
proximation can be written as 


U T {qt,q) = gls^ H T 



(3.79) 


with some function H?■ We may neglect the term q 2 — q [j from the free propa¬ 
gator, because this term does not become singular for small x. 4 After expansion 

4 However, q 2 and go would have to be taken into account when summing up the IR contri¬ 
butions to the coefficient of a 5 , since for this coefficient also less IR singular contributions are 
important. 
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of the integrand with respect to a as in Eq. (3.56) we get then integrals of the 
following type contributing to the self energy. 


ma 


5/3 


dz z 



ma 3 v n \v\ 2 n 
n — 2 


S 2 

ma°v 
n — 2 


(3.80) 


[In the last step we have used the fact that R eDx{qt, q) is an even function of t.] 
Now we see clearly that from arbitrary powers of a in the integrand we get 
contributions to the order a 3 in Re XL. The case n = 2 corresponds to the term of 
order a 3 log a, which we have evaluated already in the previous section. As we 
are interested only in contributions from the IR region, we may take oo as upper 
integration limit in Eq. (3.80), since for n > 2 we get then no contribution from 
the upper integration limit. [The cases n < 2 have been evaluated explicitly in 
Eq. (3.78).] Furthermore we see from Eq. (3.80) that from the ^-integration we 
always get a factor 
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Now we can write the function C 3 t(A) as 
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7=(V 


where b = a 3 , and the prime denotes the omission of the tree level propaga¬ 
tor. Eq. (3.82) is in fact independent of b, therefore we may simply set 6=1. 
Summing up the (Taylor) series we find 
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(3.84) 


This integral can probably not be done analytically. Numerically one finds read¬ 
ily Ji = 0.000233964448.... 

For the longitudinal part a completely analogous calculation gives 


C3L = -A(A 2 + 7T 2 ) 
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with 

1 [z 7T 2 -4 

J 2 = --^2 J dz 2 -g7-1- zm 2 Re D ' L \q^ zqo . (3.86) 

Numerically one has J 2 = 0.000250187944... 

The constant under the logarithm in the fourth line of Eq. (3.77) can be de¬ 
termined from the zero temperature limit of 03 (A), with the result 

exp < — 7 - 77 - 7 -— 97 t 6 (—5 + 6 log 2 ) 

| 12(6144 — 2567T 2 + 367r 4 — 97 t 6 ) 

+72tt 4 (-4 + 31og2) + 128 tt 2 (31 - 121og(f)) - 6144(7 - 61og(f)) 
+3456 tt 6 ( Ji + J 2 ) ) = 0.92789731088 ... (3.87) 


Let us turn now to the imaginary part. The function d 3 (A) in Eq. (3.72) can 
be determined by summing up IR enhanced contributions in a similar way as in 
the computation of 03 (A). We find 

d s(A) = ( 192 °~^ 2 5 ~ 97r4 +J 3 + [-2(Li 3 (-e- A ) + Li 3 (—e A )) - 3((3)] , 

(3.88) 
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(3.90) 


where the prime denotes the omission of the tree level propagators in the spec¬ 
tral densities. In contrast to the real part of £+, these integrals can be evaluated 
analytically rather easily, as we will demonstrate in the following. After the 
substitution z = 1 /u we can write J 3 as 


d 3 = A\ + A -2 + A 3 


(3.91) 
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, (3.94) 
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Figure 3.3: Integration contours for A 3 . 


where 0 < e <C 1. In A 3 we have used the fact the the imaginary part in the 
integrand vanishes for \u\ > 1 , which allows us to extend the outer limits of the 
integration domain from ±1 to ±00 without changing the value of the integral. 
It is easy to check that A\ vanishes for e —> 0. For A 2 we obtain 


A2 


A 3 can be written as 


1 8 48 - 7r 2 

67r 3 e 3 7r 5 e 67T 5 
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(3.95) 
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where the contour C consists of the four straight lines pinching the real axis in 
Fig. 3.3. The radius of the two small arcs in Fig. 3.3 is equal to e, whereas the 
radius of the two large arcs will be sent to infinity eventually. Since the integral 
along a closed contour vanishes if the integrand is analytic in the domain en¬ 
closed by the contour, we may replace the integral along the straight lines with 
minus the integral along the arcs. After the substitution u = re lv we expand the 
integrand for small r to obtain the contribution from the small arcs, and for large 
r to obtain the contribution from the large arcs. From the small arcs we get 
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3 67 r 3 e 3 
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and from the large arcs we get 



(3.97) 


(3.98) 
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(s') 

Figure 3.4: Re£;y at T = 0 (continuous curve), T = m (dashed curve), and T = 3 m 
(dotted curve). 


Adding up the A's we find for J 3 after taking the limit e —> 0 


J 3 


-1152 + 40vr 2 + 97r 4 
487T 5 


(3.99) 


In a completely analogous way we find for the longitudinal contribution 


12 - 7 r 2 

J 4 =-—- • (3.100) 

647r 

Inserting these results for J 3 and J 4 into Eq. (3.88), and taking the zero temper¬ 
ature limit, we find 



6144 — 2567T 2 + 367T 4 — 97t 6 7r 
864 tt 6 X 2* 


(3.101) 


as stated in the Eq. (3.77). 


3.5 Numerical results 

For intermediate energies \E — n | >m the integrals in Eqs. (3.52) and (3.62) have 
to be evaluated numerically. Figs. 3.4 and 3.5 show the real and imaginary parts 
of the self energy as a function of the energy at T = 0, T = m and T = 3 m. 
In the imaginary part we have subtracted the quantity ImS + (e = 0) which is 
divergent at finite T [see the discussion before Eq. (3.62)]. The resulting function 
ImE_)_(e) — ImE + (0) is even in e, which leads to a cusp for £ —> 0 at T = 0, while 
at finite T it vanishes quadratically for e —> 0. The real part is an odd function 
with respect to e. Again, the function is nonanalytic at £ = 0 for T = 0, but at 
finite T it is analytic at this point. At T = 0 we find 

lim ImS + (£) = —g 2 Cfm x 0.04053... (3.102) 
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Figure 3.5: Im£+(e) — ImS + (0) at T = 0 (continuous curve), T = m (dashed curve), 
and T = 3 m (dotted curve). 


R eV®\ T=0 /(g*C f m) 




(s) 

Figure 3.6: Exact expression for Re£+ | T=0 (continuous curve) compared to two per- 
turbative approximations (dashed curve: O(slogs), dotted curve: 0(s 3 loge)). 


In Figs. 3.6 and 3.7 we show a comparison of the exact (numerical) results 
at T = 0 with perturbative results as obtained from Eq. (3.77). We observe 
that the leading terms (0(elog |e|) for the real part and O (|e|) for the imaginary 
part) are already quite good approximations for e/m < 1. The higher orders 
in the perturbative expansion do not lead to a considerable improvement of the 
approximation, which indicates that the series converges well only for e/m <C 1. 

Fig. 3.8 shows a plot for the inverse group velocity [89], 



(9ReS + 

~dE 


(3.103) 


Again one can see that the leading logarithmic approximation is already quite 
close to the exact result for small s, while the hard contribution M^./(2E 2 ) 
[~ g 2 Cf/(8TT 2 ) for e <C p] dominates for larger values of s. 
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Figure 3.7: Exact expression for ImS + | T=0 (continuous curve) compared to two per¬ 
turbative approximations (dashed curve: 0 (e), dotted curve: 0 (s 3 )). 


(v- 1 - l)/(g 2 C f ) 



l°gio(e/m) 

Figure 3.8: (vf 1 — 1 )/(g 2 Cf) as a function of log W ((E — p)/m ) at zero temperature 
(from [89]). 
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3.6 Fermi surface properties 


The quasiparticle distribution function at zero temperature is given by 


n(p) 


- dpo Im- 1 -. 

7T J-oo -P0 + P + S+(po,p) - IT] 


(3.104) 


We are interested in the behavior of the distribution function for p ~ p. A possi¬ 
ble non-regular behavior is expected from the singularity of the integrand. Thus 
we consider the quantity [90] 


n(p ) 


vr L dp ° Im -p 0 + p + E+ (p 0 , p) - ip' 


(3.105) 


with 5 < //. We will neglect the imaginary part of the self energy in the follow¬ 
ing, since the imaginary part does not show a logarithmic enhancement in the 
vicinity of the Fermi surface. Thus we may write 


n(p) ^ / dpo 6(p 0 -p- Re£+(p 0 ,p)) 
J /1—8 


rv 1 
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dpo 


(3.106) 

where ca(p) is the solution of the approximate dispersion law po — p — ReS + = 0. 
We obtain 


n(p) — 


©(p-w(p))- 


1 

c>R.cS + | 

dpo \p 0 =O(p) 


(3.107) 


The second factor on the right hand side is the group velocity. In contrast to 
a normal Fermi liquid it vanishes (like l/log|ta — //|) for Co —> //. Therefore 
also the discontinuity at the Fermi surface vanishes. We remark that a similar 
behavior can also be found in condensed matter physics, in the context of high- 
T c superconductors. These systems have been termed "marginal Fermi liquids" 
[91]. 
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Chapter 4 


Anomalous specific heat at low 
temperature 

4.1 General remarks 

In this chapter we will continue our analysis of normal degenerate quark mat¬ 
ter. We will compute the specific heat, which may be relevant for cooling of 
(proto-)neutron stars at temperatures above the color superconductivity phase 
transition. For a normal Fermi liquid (see Sec. 4.2) the specific heat is linear in 
the temperature for small temperatures. We will see, however, that the specific 
heat of cold dense quark matter contains an anomalous leading term of the order 
T log T. As for the quark self energy, this logarithmic enhancement comes from 
long-range chromomagnetic fields. In this chapter we will also correct an error 
in a recent computation of the specific heat [92], which would have resulted in a 
T 3 log T behavior of the specific heat at leading order. 

4.2 Landau Fermi liquid theory 

In this section we will give a brief account of Landau's theory of Fermi liquids 
[93, 94], with particular emphasis on the specific heat. Our presentation follows 
closely Refs. [95, 96], to which we refer for further details. 

To begin with, let us consider a uniform, non-relativistic 1 gas of free spin— 
fermions. The ground state consists of a Fermi sea, corresponding to the occu¬ 
pation number 

n(k) = 2Q(k F - k), (4.1) 

where k F is the Fermi momentum. The total energy is given by 

£ = E £->(*)■ < 4 - 2 > 

_ k 

1 The generalization of Landau's theory to relativistic systems is discussed in [97]. 
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Now let us turn on interactions between the particles adiabatically. We assume 
that the quantum mechanical states of the free system are gradually transformed 
into states of the interacting system. If this is the case, the interacting system is 
called a (Landau) Fermi liquid. 

Next let us add an additional particle with momentum k (k > kp) to the free 
system, and then switch on the interactions. The eigenstate of the interacting 
system, which is obtained in this way, is called a quasiparticle. Denoting the 
energy of a quasiparticle with ek and the chemical potential with //, one finds 
that the lifetime of quasiparticles is proportional to (ek — // ) 2 [96, 98] [see also 
Eq. (3.66) above]. Therefore the concept of a quasiparticle is well-defined in the 
vicinity of the Fermi surface. Since quasiparticles are adiabatically evolved from 
fermions, their distribution function at finite temperature is given by the usual 
Fermi-Dirac distribution, 

n( k) = — - —pp -. (4.3) 

v 7 e Fk -p)/t + i v 7 

A weak perturbation of the interacting system will induce a change 6 n(k) in 
the occupation number. Landau postulated that the change in the total energy 
of the system is then given by 2 

SE = J2 4^(k) + ^ ^ k ’ k')Mk)Mk'), (4-4) 

k Z k,k' 


where is the energy of a quasiparticle, and /(k, k') is the quasiparticle inter¬ 
action. For isotropic systems the effective mass m* is defined as 


k F 

rn* 


K 

dk 


k=k ]7 


The specific heat at constant volume and per unit volume is given by 


(4.5) 




V \dT 


v 


(4.6) 


Using Eqs. (4.3)-(4.6) one finds that the specific heat of a Fermi liquid is linear 
in the temperature at low temperature (as for a free Fermi gas), C v = vT, where 
the coefficient 7 is given by [96] 

7 = ^ m*k F . (4.7) 

Landau's theory is corroborated by results of the microscopic theory, i.e. 
quantum field theory. For instance one finds also in quantum field theory that 
the specific heat of a cold fermionic system is linear in the temperature at lead¬ 
ing order in many cases. The coefficient 7 can be computed e.g. from the 2PI 

2 For simplicity we suppress spin labels. 
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effective action. For a theory with instantaneous four fermion interaction one 


finds [90] 


2tf d 

7 = -y Im X! log ( £ k - W + ReE k (w) + irj) 

k 


UJ = fl 


(4.8) 


where E k (w) is the fermion self energy. Evaluating the derivative with respect 
to u>, one finds that the coefficient is proportional to the inverse group velocity. 


(4.9) 

LO=fi, k= on shell 

It is clear, however, that Eq. (4.8) becomes meaningless if the group velocity 
vanishes at oj = //. In the previous chapter we did indeed encounter a situation 
where this happens, namely for a system with long-range interactions. In that 
case Landau's theory is not applicable, and the leading term of the specific heat 
will not be linear in the temperature. 


_i , <9ReS k 

v > 


4.3 Specific heat from the entropy 


The specific heat C v at constant volume and per unit volume can be defined as 
the logarithmic derivative of the entropy density with respect to temperature at 
constant number density. 



(4.10) 


This can be rewritten in terms of derivatives at constant T or p in the following 
way [99], 


C V = T 



(my (my 1 
yarjydfiJT 


(4.11) 


At low temperatures only the term with the entropy contributes [100], 


Cv = T (yP) +(D(T 3) - (412) 

The thermodynamic potential for QCD is given by the following functional 
of the full propagators and self energies (assuming a ghost-free gauge) [101], 

pn[D,S] = ^TV log D~ l - ^Trn D - TrlogS -1 + TrXS + $[D,S\, (4.13) 

where S is the quark propagator, E is the quark self energy, D is the gluon propa¬ 
gator, II is the gluon self energy, and 4> is a series of 2-particle-irreducible (skele¬ 
ton) diagrams. 

Using the fact that Q[D, 5] is stationary with respect to variations of D and 
S, one can derive an expression for the entropy which to two-loop order in the 
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Figure 4.1: Sum of the gluon ring diagrams. 


skeleton expansion is entirely given by propagators and self energies [24, 102], 
We may neglect the contribution from longitudinal gluons, since they are sub¬ 
ject to Debye screening and give therefore only a contribution to the normal 
Fermi liquid part of the entropy. Furthermore we may neglect antiparticle con¬ 
tributions in the fermionic sector since they would only lead to exponentially 
suppressed terms, oc exp(— fi/T). Thus we find for the entropy density (with 


P = -n/V) 



ImFEr FfeZ>r) 


ImS_|_ ReS+) + S', (4.14) 


where I) r 1 = —% + q 2 + 11/ and ,S' 4 1 = — ko + k + E + as in chapter 3. For the 
evaluation of the imaginary parts we always assume that qo and /cq are replaced 
by qo + irj and ko + irj, respectively [see Eqs. (A.9) and (A.12)]. 

In the original derivation of the anomalous specific heat in QED by Holstein 
et al. [103], only the term involving Im log Sf ] had been taken into account in the 
formula for the entropy density. While such a formula for the entropy density is 
known to be correct for standard Fermi liquid systems [90], it is not clear a priori 
how to justify it for more complicated systems, which possibly exhibit devia¬ 
tions from the Fermi liquid behavior. Therefore we will use the full expression 
(4.14) to calculate the entropy density at the order g 2 TlogT in the following. 
We will only neglect S', since it vanishes at two-loop order [102, 24], and could 
therefore only give contributions which are suppressed by an additional factor 


4.3.1 Quark part 

In Eq. (4.14) we have the following contribution from the quarks. 



x Im log (—ko + k + S + ) — ImE + Re 


1 


ko + k + 
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Im log (—ko + k) + ReS + Im 


1 


-feo + k 


(4.15) 


where we have performed an expansion with respect to £+, keeping only the 
free term and the term corresponding to a single quark self energy insertion. 
Diagrammatically, the latter part corresponds to the sum of the gluon ring di¬ 
agrams, see Fig. 4.1. The free term gives the particle contribution to the free 
fermionic entropy density, 

S^~NN f ^l. (4.16) 

In the last term in Eq. (4.15) the factor Iml /(—ko + k ) forces the self energy to 
be on the light cone. From Eq. (3.67) we obtain the following approximation of 
Re£ + at leading order 


ReX+ ~ 


r£l 

127T 2 


Oo - n) 


log 


M 


\k 0 - n\ 


+ h 


f k 0 - fi Y 
V T )\ ’ 


(4.17) 


with M cx f/efr/i, and h is some function. From the Fermi-Dirac distribution in 
Eq. (4.15) we see that ko — // is of the order T. At leading logarithmic order we 
may therefore drop h in Eq. (4.17), since it gives only a 0(1) contribution to the 
result. This contribution would fix the scale under the logarithm, but this scale 
will be determined anyway in Sec. 4.5 using a slightly different approach. Thus 
we get 


C NLO 

‘Aq) 




(4.18) 


We make the substitution k = Tz + fj,. The integral is dominated by small values 
of z, therefore we may send the lower integration limit to —oo. Then we obtain 
at order T log T 


cNLO 

Wq) 


9cS N 9d 2T 

367r 2 


log 



(4.19) 


This result agrees with the one of Holstein et al. [103] after correcting a factor of 
4 therein, as done previously in [104], 


4.3.2 Gluon part 

The gluon part <S( g ) is given by the first line of Eq. (4.14). It can again be inter¬ 
preted in terms of ring diagrams. Fig. 4.1. Using the relation [24] 


i ( ImII 

Im log D = arctan ————r 
yRei9 -i 


■ ne(qo)9 (-R eD x ) , 


(4.20) 
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we write S (g) = 5 (cut) + 5 (n) + «S( po i e )/ with 


5(cut) — 2 N g 


s 


£(n) - -2 N g 


^(pole) - 2 Ng 


d 4 Q dn b (q 0 ) 
(2vr) 4 dT 
d 4 Q dn b (q 0 ) 
(2vr) 4 8 T 
d 4 Q dn b (q 0 ) 
(2vr) 4 8 T 


arctan 


[ml [7 


<?o - q 2 ~ ReII T 


IiiiIIt’ Re 


% — T ~ Hr 


(4.21) 

(4.22) 


vre(g o )0 (<7o ~ 9 2 ~ Mr) • (4.23) 


For the cut term we use the approximation go ■C q, because it can be checked 
that including terms of higher order in go would only produce terms of higher 
order than TlogT (see Sec. 4.5). In this region we have from Eqs. (3.29) and 
(3.31) 


Fir — —i 


_2 ,, 2 „ 
■9el tM 9o 


Anq 


Introducing an UV-cutoff q c for the moment, we obtain 


S(cu„ = 2| / rf® 


ric 


dq q‘ 


9o 


dn b (q 0 ) 

dT 


arctan 


ffeffAo 

47rg 3 


(4.24) 


(4.25) 


In order to evaluate this integral we make the substitution y = qo/T, x = 
4tt q 3 /(g 2 s y 2 T). Keeping only the term of order TlogT, we get the (finite) re¬ 
sult 


C ^9 2 eS NgH 2 T (M’ 
5( C ut) - 36vr 2 lQ g ( T 


(4.26) 


The determination of the constant M' would require a more accurate calculation 
including the longitudinal mode (see Sec. 4.5), which would show that M' oc 

Next let us evaluate 5 (n )- Following similar steps as in the computation of 
S(cut), we find 


5 ( n) - 2Ng 


dqp f qc dq q 2 dn b (q 0 ) ( g 2 s g 2 qo 
2vr 


'90 


9lsN g g 2 T 

367T 2 


log 


2vr 2 dT 


NT 
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Airq 
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9eft A' 2 90 N 2 
47T q 


(4.27) 


We observe that at order T log T this expression just cancels the contribution 
from Eq. (4.26). 

Finally let us consider the pole part. In the HDL approximation we have at 


low temperature ~ 211 t, and therefore 


dS, 


T- 


(pole) 


dg 


dfi 


= -2 N„ 


= -4v tN„ 


d 4 Q dn b (q 0 ) 
(2vr) 4 dT 
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(% ~ 9 Mg o) R9 0 -r- ReFR , (4.28) 
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Figure 4.2: Illustration of the inequality (4.31), with the plasma frequency w p normal¬ 
ized to 1. Continuous line: c ot(q), dashed line: Q(lo p — q) + q@(q — w p ), dash-dotted 
line: oj p + q. 


where we have discarded contributions ~ T 3 which are negligible in the low- 
temperature limit. Using [76] 


e(q 0 ) 6 [ReDj . 1 ) = Z T (q) (%o - u T (q)) ~ %o + w T (q))) (4.29) 


we find 


h 


(pole) 2 IV f j 

dp IT 2 


dq q 2 dnh ^T(q)) ( a>T ( g )2 _ q 2 )z T (q). (4.30) 


We would like to give an estimate for this integral. With the plasma frequency 




= 9 eS P / (tt \/3) we have the inequalities 


(see Fig. 4.2), and 


Q(u p - q) + qQ(q - u p ) < u T {q) <w p + q 


1 


z T{q) < tt 
2 q 


(4.31) 

(4.32) 


(see Fig. B.2). Since we assume T <^u p o c g e sp, we can therefore estimate 
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(4.33) 


and 


dq q ' 


; dn b (u) T (q)) 
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dqq(u p + q) 2 dnb ^ + ^ 


dqq(u p + q) 


2 9n b {q) 


}H e -„ p/ T 

24 T 2 


(4.34) 


Apart from terms ~ T 3 which are dropped in the derivative with respect to 
//, this crude estimate (which we shall refine in Sec. 4.5.5) shows that the pole 
contribution is exponentially suppressed (essentially because of lot > cc p ). 


4.3.3 Result 

From the previous subsections we find the following result for the entropy den¬ 
sity at low temperature. 


cS = 5 (g) + 5 (q) ~ NN f ^- + 9cS ^ 2 T log (^ 


(4.35) 


From Eq. (4.26) we see that the T log T term can also be obtained by starting 
only with 

n at f d4 Q n -i /A ~,, 


S~NN f t A _ 21 V 9 J ' Imlog D f (436) 

This formula corresponds to integrating out the fermions, as has indeed been 
done in the approach of [100, 105]. On the other hand, we see from Eqs. (4.15) 
and (4.19) that one also gets the correct result by using the purely fermionic 
expression 


-Imlog Dj, 1 . 


S ~ -4 NN f J — Imlog(-fc 0 + k + ReS + ), 


(4.37) 


which justifies the starting point of [103]. 


4.3.4 A note on the Sommerfeld expansion 

Since the inverse group velocity diverges on the Fermi surface, one might won¬ 
der whether the coefficient of the linear power of temperature in the specific 
heat should also diverge. Of course we have seen from explicit calculations that 
this is not the case. In order to see in which way a straightforward application 
of Fermi liquid techniques goes wrong in our case, let us consider again the 
entropy density. From Eq. (4.15) we have 

S NLO = -\NN } dk 0 k f nf{j lr /X) ReF+| T =o. (4.38) 

Jo dT 

A naive Sommerfeld expansion, assuming (erroneously) that ReE+|T=o is ana¬ 
lytic at ko = //, would lead to 

sNLO - - NN fim; « ReS +> = ” sVT ' < 439 > 
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Then we would have for the part of the entropy density which is linear in tem¬ 
perature 

Sun. ~ NN f ^( 1 + Voo”), (4.40) 

Of course this equation is wrong, since the assumption that ReS + |r=o is analytic 
at ko = n is incorrect. The explicit calculation (see above) shows that instead of 
a linear term with a singular coefficient, one obtains a T log T term with a finite 
coefficient. 


4.4 Specific heat from the energy density 


The energy density can be obtained from the expectation value of the energy 
momentum tensor, 

U = y J d 3 x(T 00 (x)), (4.41) 

and the specific heat is then given by 


C 


v — 



(4.42) 


Here the temperature derivative has to be taken at constant particle number 
density, in contrast to the calculation of the low temperature specific heat in the 
previous section, where all temperature derivatives had to be taken at constant 
chemical potential, see Eqs. (4.12), (4.14). Of course, both methods will ulti¬ 
mately yield the same result. 

In [92] the specific heat is computed using the following formula for the total 
energy density. 


U = 2 


dko 


d 3 k 

(27r) 3 


rif(ko - 


fi)k 0 p + (k 0 ,k), 


(4.43) 


where p + is the spectral density of the positive energy component of the quark 
propagator (see below). It should be noted that this formula is incorrect even for 
a theory with only instantaneous interactions of the type 


Hint = 7 } J d 3 x d 3 x' v4(x, t)V aa ',pj3'{ x - x / )^ j g/(x / , t)Vv( x , i), (4-44) 

in which case the correct formula reads [106] 


U = 2 


dko 


d 3 k 

(27t) 3 


nf(k 0 - 


k) 2^0 + k)p + (k Q ,k). 


(4.45) 


The anomalous behavior of the specific heat comes from dynamically screened 
interactions, whose non-instantaneous character cannot be neglected. It would 
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be rather difficult to generalize Eq. (4.45) directly for non-instantaneous inter¬ 
actions, because one would have to use an effective Hamiltonian which is non¬ 
local in time. If one naively tries to compute the anomalous specific heat from 
Eq. (4.45), the coefficient of the logarithmic term turns out to be wrong by a fac¬ 
tor . Therefore, we will follow a slightly different approach, using the energy 
momentum tensor of QCD without integrating out the gluons. 

The energy momentum tensor can be written as a sum of three pieces. 


' J (q) + 4g) +J W’ 


(4.46) 


corresponding to the quark part, the gluon part, and the interaction part. The 
contributions of these parts will be evaluated in the following subsections. We 
will neglect gluon self interactions and ghost contributions, since they give only 
higher order corrections at low temperature. 


4.4.1 Quark part 

The quark part is given by 




(4.47) 


where we have written explicitly the sum over flavor space. This is the (only) 
contribution which is taken into account by Boyanovsky and de Vega [92], Let 
us repeat their calculation, but for simplicity without the renormalization group 
improvement of the quark propagator proposed in [107]. As for the entropy 
density is is sufficient to take into account only the positive energy component 
of the quark propagator. Thus we find 


q) = 2NNf / dko 


d 3 k 

( 2 ^ 


nf(k 0 - p)k 0 p+{k 0 ,k), 


(4.48) 


where the spectral density is defined as p + = ilmS , +- 

In order to obtain the specific heat from Eq. (4.42), we have to determine first 
the temperature dependence of the chemical potential from the condition 


dM _ dM dp dM 
~dT = ~df + ~dT~dp = ’ 

where the particle number density M is given by 

d 3 k 


(4.49) 


/ [■ ppU 

dko J J^ 3 n f(ko - d)p+(k 0 , k) (4.50) 


(up to anti-particle contributions). We expand J\f with respect to g, 

J\f = A/q + 2 T • • • 


(4.51) 
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The free contribution A/o is given by 

A/" 0 = NN f 


(4.52) 


In A /2 we are only interested in contributions which contain log(M/T). Such 
terms arise from infrared singularities from the gluon propagator, which are 
dynamically screened. This corresponds to scattering processes of quarks which 
are close to the Fermi surface. Therefore the anomalous terms come from the 
region k ~ ko ~ p, where X + is given by (4.17) at leading order. At leading 
logarithmic order we may again neglect the function h in Eq. (4.17). Subtracting 
the temperature independent part, we find then 

/ r (fib 

dk 0 / 3 (n f (k 0 - p) - 0(p - k 0 )) 5(k 0 -k-ReT, + ) , 

J ^> 0(9 e 2 ff) 

(4.53) 

where we have approximated the spectral density with a delta function, since 
the imaginary part of X + is negligible compared to its real part. The integration 
can be performed easily, with the result 


2 K r 9eS N gfi T2 , 

g eS ^2 - C 3 6 ^ 2 log 


(4.54) 


We notice that this result is consistent with the result for the entropy density, Eq. 
(4.35). Now we can solve Eq. (4.49) at low temperature, 

dT 3/j 18 & \T J v ’ 

The approximate solution to this differential equation is given by 

MT) = «(0)(l-^5-f^log(f)). (4.56) 

Eqs. (4.55) and (4.56) correctly reproduce the beginning of the perturbative ex¬ 
pansions of Eqs. (2.37) and (2.38) in Ref. [92], 

For the specific heat we obtain from Eq. (4.48), following the same steps as 
in the calculation of dJ\f / dT, 


l4 T ) - l-t(O) 1 - 




(4.55) 


(4.56) 


C„„ = NNj/i 2 T + NN,^ + l log 


(4.57) 


Using Eqs. (4.55) and (4.56) we find that the T log T-terms cancel, 

C v{q) ^ NN f ^f, (4.58) 

as stated in [92], We should emphasize that this cancellation has nothing to do 
with the nonperturbative renormalization group method which is employed in 
[92], It should also be noted that this renormalization group method has been 
criticized recently [108], because the neglect of the /3-function in the renormal¬ 
ization group equation used in [92,107] does not seem to be justified. 
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4.4.2 Gluon part 

We now turn to the gluon part of the energy density, which has been explicitly 
neglected in Ref. [92], This is given by 

Tg = i (E“ • E a + B“ • B a ) . (4.59) 

Neglecting gluon self interactions, and keeping only the transverse part of the 
gluon propagator, we obtain 

W( g ) - 2 N g J j^ytMqo) Im ((c/o + Q 2 ) D t ) • (4.60) 

The pole contribution to this integral is again exponentially suppressed, there¬ 
fore we only have to consider the cut contribution. At low temperature we may 
neglect the temperature dependence of the gluon self energy. Subtracting a less 
infrared sensitive contribution not involving the Bose distribution (correspond¬ 
ing to the non-rife-part in Sec. 4.5), we find 

,4 n 

(n b (q 0 ) + O(-q 0 )) («o + « 2 )- 2 2 .r ( 4 - 61 ) 

W + (%!») 

After the substitution qo = Ty, q ' = <jf f( fi 2 Tx/(47T) the integral can be readily 
done, with the result 

W (g) ^ log (M/T), (4.62) 

which gives the following contribution to the specific heat at order T log T, 

<463) 

Again the determination of the constant M would require a more accurate cal¬ 
culation, similar to the one in Sec. 4.5. 

4.4.3 Interaction part 

The interaction part is given by 

i?L)=s'£^P 0A ° T ‘*r < 464 > 

/ 

The expectation value of this term is essentially given by the resummed gluon 
ring diagram as in Fig. 4.1. However, here only the longitudinal component 
of the gluon propagator appears in the loop. This mode is subject to Debye 
screening, so it can contribute only to the normal Fermi liquid part of the specific 
heat. 
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4.4.4 Result 


From the previous subsections we see that the only contribution to the specific 
heat at order g 2 T log T when calculated along the lines of Ref. [92] comes from 
the gluon part, Eq. (4.63). While this confirms the observation of Ref. [92] that 
the quark contribution of order g 2 T log T cancels against a similar term from the 
temperature dependence of the chemical potential at fixed number density, it 
shows that the neglect of the gluon contribution to the energy density in [92] is 
not justified. 

Our result for the specific heat from the energy density is consistent with 
the result for the entropy density in Sec. 4.4. For the entropy density we found 
at leading logarithmic order that it was sufficient to take into account only the 
contribution of fermionic quasiparticles. For the energy density, on the other 
hand, the purely fermionic formula (4.48) is not sufficient to give the correct 
result. 


4.5 Higher orders in the specific heat 


In this section we will evaluate higher terms in the low temperature specific heat 
which go beyond the leading logarithmic approximation. A convenient starting 
point is the following expression for the pressure, which becomes exact in the 
limit of large flavor number Nf [109,110, 111], 


P = NN 
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1 D -1 
--Imlog-^ 

2 q 2 ~ qo 


7 ;Irn log D v A 
+ 0( 5 V),( 4.65) 


where the inverse gauge boson propagators are given by D T 1 = q 2 — + H r + 
n vac , Dl l = q 2 - ql + n L + n vac , and D~£. = q 2 - q 2 + n vac . The first line of Eq. 
(4.65) is simply the pressure of a free fermionic quantum gas, and the remaining 
part is essentially given by the sum of gluon ring diagrams in Fig. 4.1. At finite 
N f , Eq. (4.65) with II including also the leading contributions from gluon loops 
still collects all infrared-sensitive contributions up to and including three-loop 
order [112, 113]. We shall however find that all contributions from gluon loops 
to II are suppressed with (T///) 2 in the specific heat, and will thus be negligible 
for T <C ^ compared to the terms we shall keep. 

In the following we will always drop temperature independent terms in the 
pressure, since they do not contribute to the specific heat at low temperature. We 
will refer to terms involving the bosonic distribution function rib as "n^-parts", 
and to the remaining ones as "non-n/,-parts". 
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From the explicit form of the self energy (see Sec. 3.2 and [77]) we see that 
the (/-integration naturally splits into three regions: q < go, qo < q < 2y — go, and 
2/i —(/o < q < 2p + g 0 . We will consider these regions separately in the following. 


4.5.1 Transverse contribution, region II 

The n/,-part of the contribution of the transverse gluons to the pressure is given 
by 


Pt 0 f d 3 q [°° dq 0 2 2 , tt _i_ tt \ 

TT = ~ 2 7TTT / -n 6 (g 0 )Imlog(g -q 0 + U T + n vac ). 

J (2vr)' 5 J o n 


The cut contribution from region II is given by 
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nft(go) arctan 


Imllr 


g L ~~ 9o + heffr + ReII v 


(4.66) 


(4.67) 


As long as T C/i it is sufficient to take the self energy at zero temperature. 

The Bose-Einstein factor and the leading term in the gluon self energy set the 
characteristic scales, 

90 ~ T, q ~ (<&/i 2 T) 1/3 . (4.68) 

We want to perform an expansion with respect to a parameter b defined as 


b m ,/s 

\ .9eff /I / 


(4.69) 


It turns out that the following approximation of the gluon self energy is sufficient 
through order T 3 log T in the entropy density. 


ReU T (q 0 ,q) ~ 

7T Z 

Imn T (g 0 ,g) ^ 


q 2 3 q 4 J ’ 


q q 3 



(4.70) 

(4.71) 


The first two terms in both lines are the leading terms of an expansion of the 
HDL self energy in powers of go- Naively counting powers of b in the integrand 
one would conclude that only these terms are responsible for the coefficients 
of the expansion in b, while the other terms in the self energy should give only 
terms which are suppressed with additional powers of g in the pressure. In prin¬ 
ciple this is correct, but it turns out that it is necessary to keep also the last term 
in Eq. (4.71), which is beyond the HDL approximation. As we well see in Sec. 
4.5.3 this term gives a contribution in the hard region, where the approximation 
q <C n is not justified. 

Introducing dimensionless integration variables x and y via 

qo = b 3 g e s^y, q = bg e sn(x/( 4vr)) 1/3 , (4.72) 
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we find after expanding the integrand with respect to b, 
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(4.73) 


3ir 2 x(x 2 + y 2 ) 3 

with .x„un = 47r& 6 y 3 and x max = 47r(2— b 3 g e sy) 3 / (bg e g) 3 . In the coefficients of this 
expansion we have written down only those terms which do not ultimately lead 
to terms that are suppressed by explicit positive powers of y e ff • The integrations 
are now straightforward, and we find 3 
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(4.74) 


The evaluation of the constant ct is a bit more involved because one has to sum 
up an infinite series of contributions from the infrared region. This will be done 
in Sec. 4.5.7. 


4.5.2 Longitudinal contribution, region II 

The rtft-part of the contribution of the longitudinal gluons to the pressure is given 
by 
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9 J ( 27r ) 3 Jo tt 
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n 6 (y 0 )Imlog 1 


q 2 -qo 


(4.75) 


As in the previous section the ^-integration decomposes into three parts. In the 
second region (qo < q < 2y — qo) the characteristic scales are now 


qo 


q ~ g e sg- 


(4.76) 


In a similar way as in the previous section, the gluon self energy can be approx¬ 
imated as 
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(4.78) 


3 For the y-integration one needs the following formulae valid for a > 1: 

.a — 1 poo „.oc — 1 i 

dy- 

/ o eV - 


POO OL — 1 pOO OL — 1 1 

J y = r («K(«), j dy y eV _l V = r'(a)C(a) + r(a)C '(a). 
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We introduce dimensionless integration variables y and z via 

qo = b 3 g e gyy, q = g e gyz/n. (4.79) 


Then we find after expanding the integrand with respect to b, 


N n 


4 4 

ffcffV 

16ir 2 


dy 


ev -1 


/: 


dz 


min 

4^ 


6 yz(-4ir 2 (l+z 2 )+g 2 s z 4 ) 12 y 3 (ir 2 - 12(z 2 + 1)) 


vr 4 (l+z 2 ) 2 




3z(l 


>2\3 


+ 0(b 


18 > 


(4.80) 


with z m i n = b 3 yir and z max = (2 /g e g — b 3 y)ir. In the coefficients of this expansion 
we have written down only those terms which do not ultimately lead to terms 
that are suppressed by explicit positive powers of g e g. The integrations are now 
straightforward, and we find 


pii 

1 t, 4 4 

JT g * 


b 6 


+ 


487r 2 

7T 2 (12 — 7T 2 ) 
240 


i + bglgl 


b 12 log b + CLb iZ + (9(6 i8 log b) 


t.12 


18 i 


(4.81) 


As in the transverse sector, the determination of the constant cl requires sum¬ 
mation of contributions from the infrared region, which will be done in Sec. 


4.5.7. 


4.5.3 Contributions from the upper integration limit 


In this subsection we will show that the HDL approximation for the gluon self 
energy is sufficient to produce the series in b, apart from the order //’, where 
precisely two additional terms in the gluon the self energy are necessary, one in 
the transverse part and one in the longitudinal part. 

First let us focus on the transverse part. It is convenient to multiply the 
denominator and the numerator of the argument of the arctangent in Eq. (4.67) 
with a factor x 1 ' /3 / (y 2 fr /.x 2 /r). Dropping numerical factors, and also powers of y, 
we have then the following effective power counting rules: 


2 2 
Seflfh 



(4.82) 

(4.83) 

(4.84) 

(4.85) 


where n, r, s.t G No- Here the third line corresponds to a HDL contribution, and 
the last line corresponds to terms beyond the HDL approximation. At first sight. 
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we would expect that contributions from the last line (when r + s + t > 0 ) could 
only produce terms which are suppressed with powers of g e s in the pressure. A 
more careful analysis is required, however, since the upper integration limit is 
proportional to 1 /gfg, which might invalidate the naive power counting in the 
integrand. 

After the Taylor expansion with respect to b, the term of order b 6+2M in the 
integrand will have the following general structure for M > 1, 


4 03 


4 4 71 

9eS9 0 


2 M 


( X 2 + y 2 ) M ’ 


(4.86) 


where P 2 m is of the order b 2M . We note that P 2 m is a polynomial of degree 2 M 
in the expressions on the right hand sides of Eqs. (4.82)-(4.85). This can be seen 
from expanding first ReFT/ and [ml fy, and then the arctangent, with respect to 
b. Each term in this polynomial contains at least one factor which comes from 
ImTlr, and at least one factor which comes from q 2 — cfc + RcTTr- Explicitly, we 
find contributions to the pressure of the following type from the region of large 
x in (4.86 ), 4 


U M ■= ffeSM 4 


'•l/toeff ) 3 ^+2 M 


dx- 


,2M 


a+ — 
-X 3 


n 

i= 1 


( 9 , 


2 e^) di 


n 

3 = 1 


(9efiX) fj , 


(4.87) 


where x a+ 3 comes from the free gluon propagator, the first product comes from 
Renr, and the second one from Iml R. The parameters a, 6, 7 , 5, Ci, d u e ; , and 
fj are natural numbers (including zero). The properties of Pom listed above lead 
to the following set of restrictions for these parameters. 


7 <5 

4/3 + 2 ^2 Ci + 2 X] e 3 = 2 M, (4.88) 

i =1 j= 1 

o + /? + 7 + <5 = 2 M, (4.89) 

a + (3 + 7 > 1 , (4.90) 

8 > 1. (4.91) 

We can eliminate a and (3 in Eq. (4.87) with the help of Eqs. (4.88) and (4.89). 
After performing the integration we obtain 

Um ~ b v g%, (4.92) 

where the exponents are given by 

7 5 

v = 3(l + M + ' r + 8-^d i -J2fi)’ ( 493 ) 

_ i= 1 3 = 1 


4 The upper integration limit in Eq. (4.87) corresponds to q m ax = 2/x. Actually there exists no 
sharp upper cutoff if one takes the gluon self energy at finite temperature. There is only a "fuzzy" 
cutoff which is smeared over a region of size Ag ~ T. However, this effect gives only corrections 
of order T //j ~ g e ab 3 . 
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(4.94) 


7 <5 

w = -3 + M + 37 + 3(5 - ^2 d i Y,fr 

i =1 1=1 

We are interested in contributions with w < 0. Therefore we obtain from (4.91) 
and (4.94) 

, 7 S 

3 (3 — M -3 7 + E + E /i) > 5 > !• (4-95) 

*=i i=i 

Furthermore we get from (3 > 0 using Eq. (4.88) 

7 <5 

M - > 0. (4.96) 

i=l 1=1 

Combining the inequalities (4.95) and (4.96) we get 

7 5 

3 7 + ~ <k) + “ /?') - °- ( 4 - 97 ) 

*=i l=i 

Comparing (4.85) and (4.87) we find c t — d % > — 1, and therefore 

7 

- di) > -7- (4-98) 

2=1 


In order to proceed we have to take a closer look at the structure of Imllr. In the 
explicit zero temperature expression 


Imfl^ 


9eS ( 12 Ao 0 12p 2 g;j 2 q* ql 

487T \- —~ 3qqo + ^- + T + 7 


(4.99) 


we have e-j = fj = 0 for the first term, ej = fj = 1 for the second one, and 
e j — fj > 0 for the three remaining ones. Finite temperature corrections to 1 11 i IT / 
also have e 3 — fj > 0. Because of (4.98) and 7 > 0 we conclude that the inequality 
(4.97) can only be satisfied for 7 = 0 and for contributions which come from the 
first two terms of Iml 1 / in Eq. (4.99). We write 


S = Si + 62 , 


(4.100) 


where 61 and 82 denote the numbers of factors in Eq. (4.87) which come from the 
first and the second term of Imllr, respectively. We note that 

S 8 

Y, e i=Y,fi = h- < 4 ' 101 > 

1=1 1=1 

Therefore the inequalities (4.95) and (4.96) are reduced to 

\(3-M + 5 2 )>5 1 + 8 2 >l, M - 82 > 0. (4.102) 

O 
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The condition a > 0 gives with the help of Eqs. (4.88), (4.89) and (4.101) 

3 M - 2(5, - S 2 > 0. (4.103) 

The last inequality to be satisfied is (4.90), which gives (using Eq. (4.89)) 

2 M -Si- S 2 > 1. (4.104) 

The only solution of the system (4.102)-(4.104) with non-negative <5, and S 2 is 

M = 5 2 = 1, Si = 0. (4.105) 

We conclude that the only contribution from the upper integration limit that is 
not suppressed with powers of g e g arises from the following set of parameters 
in Eq. (4.87), 


(M, a, (3, 7 , S) 


(1,1,0,0,1), (ei,/i) = (l,l), 


(4.106) 


which gives v = 6 and w = 0 for the exponents in Eq. (4.92). Explicitly, we have 
the following contribution to the pressure in this case. 


4 4 

9e 


s 


1 /(b/eff) 3 


dx b 8 g 2 eS x~ 1/3 


9e ft I 1 


4,6 


(4.107) 


We have thus shown that this is the only contribution in the whole series in b 
(for M > 1) which comes from the upper integration limit and which is not 
suppressed with powers of g e g. We conclude that the naive power counting, as 
inferred from (4.82)-(4.85), is almost correct, which means that HDL self energy 
is sufficient to produce the series in b, with the sole exception of the term of order 
if’. At this order we get a contribution from a term in I 1111 1/ which is beyond the 
HDL approximation [the second term in Eq. (4.99)], as shown in Eq. (4.107). 

The conclusions of this section can also be understood in a more direct way. 
With the definition II). ; := IT / — II 7 - we have 

Imng } ~ q > q 0 , (4.108) 

l07T 

Imn® ^ 0, q -> q () . (4.109) 


The expression in the first line is precisely the second term in Eq. (4.99). It gives 
the following hard (two-loop) contribution to the pressure, 

1 

1 r°° ImTT^ a 2 „ 2 T 2 

= ~^l "^ 8 ^' (4110) 


pII,non—HDL pII,non—HDL 

1 (J 1 


T =0 
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This is precisely the contribution which comes from the term proportional to 
g 2 ff in Eq. (4.73). Since this two-loop expression is sensitive to hard scales, it 
is not surprising that the HDL approximation is not sufficient in this case. On 
the other hand, we see from the power counting rules (4.82)-(4.85) that the HDL 
approximation is certainly good enough to produce all the terms in the series in 
b that come from soft scales. Of course one has to resum diagrams at soft scales, 
as we have done it in Sec. 4.5.1. 

The arguments of the preceding paragraph can also be applied to the longi¬ 
tudinal case. As in the transverse case the only contribution that is sensitive to 
hard momenta appears at the two-loop level. Again we define 1 Y) ; := IR — IR, 
which gives 

Imnf ~ - ge g^° , q > go, (4.111) 

Iml4 2) -> 0, q -> g 0 - (4.112) 


We have the following two-loop contribution from hard momenta. 


N, 


a l 


pII,non—HDL pII,non—HDL 

r L r L 


= ~2t^J 0 dq ° nb ^ 


T =0 

2[i 


2 Tmn^ 

dqq 


2rp2 


Qo 


T 




geffM T 

48vr 2 


(4.113) 


All the other terms in the series in b are generated by the HDL approximation. 


4.5.4 Cut contributions from regions I and III 


In the first region the imaginary parts of the complete self energies (ER + Il vac 
and IR + II vac ) vanish, therefore there is no cut contribution to the specific heat. 

In the third region we have go ~ T, q ~ 2fi and dq ~ go ~ T. If we take 
the gluon self energy at zero temperature we can therefore make the following 
estimates for the transverse part. 


pill i roo 

T =--J dq 0 


r-2/i+go 


N„ 


7T 3 Jo 


p 2 T 2 arctan 


I2fi-q 0 


dq g 2 n;,(go) arctan ■ 


•9eV /T 




4 4 7 9 

9eSH X 9eS0 


ImII r (g 0 , g) 

ql + Ren T (g 0 ,g) 

(4.114) 


and for the longitudinal part. 


P m o^T 2 

—— ~ g?T 2 arctan — 

N g ** p 2 


geffP 4 X g 2 s b 12 . 


(4.115) 


These are contributions which are suppressed by explicit powers of g e g. If one 
takes the gluon self energy at finite temperature, the sharp boundaries of region 
III will be smeared by an amount of order A q ~ T, but this does not modify the 
above estimates. 
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4.5.5 Pole part 

The pole contribution comes from region I, i.e. q < qo- In this region the HDL 
propagator has single poles for qo > u p at real qo = u.'rjJ.q), which allows us to 
carry out the go-integration. For the n/,-part we find 



/ 1 _ e —UT(q)/T \ 

\ l-e-i/ T ) 


(4.116) 



We observe that all pole contributions are exponentially suppressed for T <C lo p , 
with the only exception of a T 4 contribution from the term Im log(g 2 — gg) in the 
longitudinal sector, which contributes the equivalent of an ideal-gas pressure of 
one bosonic degree of freedom, but with negative sign. 
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4.5.6 Non -n b contribution 

Integrals without n b are less IR singular than the corresponding integrals with 
rib, since the Bose distribution behaves like T / q 0 for small q$. This means that 
within our perturbative accuracy no resummation is necessary for the non-n;, 
terms, and it is sufficient to consider the strictly perturbative two-loop expres¬ 
sion for S non _ nb . We can determine this contribution by the observation that at 
two-loop order also the rib part is IR safe and given by 


p2—loop 
*n h 


N n 


d 3 q 

(27r) 3 



—n b (<?o)Im 
7 r 


2IIt + TR \ _ 
Q 2 ~Qo ) 


d 3 q 

(27r) 3 


dq 0 


IT 


nb(qo) 


l ^5{q-q 0 )ReU G + V^- 


—-y Iml b; 

% 


, (4.122) 


where Ik; = 2I\j + 11/,. In order to evaluate Eq. (4.122) we need Rene on the 
light cone. From Eq. (3.10) we find 


ReII/;|q 0= g 




dk k n(k) = g 2 s 



(4.123) 


Therefore we get the following contribution from the real part of the gluon self 
energy. 


p2—loop, Re 
*n b 


N n 


n 2 ii 2 T 2 

9 effb 1 
247r 2 


+ 0(T a ). 


(4.124) 


Next we want to evaluate the contribution with Imn^ in Eq. (4.122). First we ap¬ 
proximate the gluon self energy with its zero zero temperature version. (We will 
check the validity of this approximation below.) The main contribution comes 
from region II, where 


-o—o Im n G 

r - do 


9cff9o 
47 rq 


(4.125) 


Inserting this expression into Eq. (4.122) we find 


p2—loop, Im 
*rib 


N n 


9> 2 T 2 

24t r 2 


(4.126) 


At this point we are able to check explicitly that at this order including finite 
temperature corrections into the gluon self energy does not change the result. 
From the integral representation of IrnEk; [Eqs. (3.13)-(3.16)] we find 


P, 


2—loop, Im 2 poo roo 

- = -|^ J o dq 0 I dqqn b {q 0 ) 


JO 

, 7 "70 , 10 + 1 

Q(q — q 0 ) / dkn(k) — @(qo — q) / dkn(k ) 

/ 7+70 / 70-7 

J ~2~ 


(4.127) 
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(Here we have dropped the vacuum contribution in Iml \q, since it would give 
no contribution in the following.) After doing the ^-integration we find for the 
second order susceptibility 


1 d 2 H,?“ loop ’ Im 


Nn 


dp 2 


^=o 


ffcff 

4vr 4 r 



dqqn b (q 0 ) 


sinh (-^p) sinh (§f) 
[cosh (-§p) + cosh (^)] 2 


(4.128) 


We substitute q = 2 T log a and qo = 2 T log 3, which gives 

1 d 2 P 2 - loop ’ Im _2 gl s T 2 r™ r°° (P 2 - 1) log/3 _g 2 eS T 2 

N g dp 2 ^ vr 4 .A Q .A P (a + 3) 2 (l + a3) 2 127T 2 

(4.129) 

Comparing Eqs. (4.126) and (4.129) we see that at this order the gluon self energy 
at zero temperature is sufficient to give the correct result. This could have been 
expected, of course, since finite temperature corrections to the gluon self energy 
are suppressed with (T/p) 2 . 

Adding up the two contributions (4.124) and (4.126), we find that they cancel 
precisely at the order g 2 B p 2 T 2 , which means that 

p2- lo °p _ Q (4.130) 


Therefore the two-loop non-n^, contribution is equal to the standard perturbative 
result at order g 2 B fj, 2 T 2 [114,115], 


P, 


non— 


N n 


9c IqJ + ° 0?eff/A S 2 eS T 4 ) ■ 


(4.131) 


4.5.7 Evaluation of the coefficient of T 4 


In this section we want to compute the complete coefficient of T 4 in the pressure. 
To this end we will perform a systematic summation of IR enhanced contribu¬ 
tions, in a similar way as in Sec. 3.4.2. 

First let us consider the transverse part. From the terms which are explicitly 
shown in Eq. (4.73) we find the following contribution to the coefficient of h 12 in 
the pressure. 


-(i) 


1 


C-rp - 


810vr 2 L 


- 36vr 2 + 


(vr 2 - 8) ( 248 - 967 ^ - 9tt 2 + 48 log(4yr) + 8640 


C'(4) 


7T 


(4.132) 


Evaluating explicitly the next few terms in the expansion of the integrand, one 
finds additional contributions of order 6 12 . They arise from the fact that the x- 
integrations would be IR divergent, were it not for the cutoff x m i n oc //’. There¬ 
fore also terms which are formally of higher order in the integrand contribute to 
the order 6 12 in the pressure. 
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Since x m i n depends on g e g only through b, and since we are not interested 
in terms in Pf which contain g e g explicitly, it is sufficient to take the HDL self 
energy in the following. With the substitution (4.72) the transverse gluon self 
energy can be written as 

II T = g 2 eS y 2 H T , (4.133) 


with some function Hr■ In this section we may neglect the term q 2 — q 2 } from 
the free propagator in Eq. (4.67), because this term does not become singular for 
small x. After expansion of the integrand with respect to b we get then integrals 
of the type 



b 2 y\ n b 12 y 3 
x 1 / 3 ) n — 3 


(4.134) 


Now we see clearly that from arbitrary powers of b in the integrand we get 
contributions to the order h 12 in Pj. The case n = 3 corresponds to the term 
of order b 12 log b, which we have evaluated already above. As we are interested 
only in contributions from the IR region, we may take oo as upper integration 
limit in Eq. (4.134), since for n > 3 we get then no contribution from the upper 
integration limit. [The cases n < 3 have been evaluated explicitly in Eq. (4.132).] 
Furthermore we see from Eq. (4.134) that from the ^-integration we always get 
a factor 


dyy 3 

e''J - 1 


7r 

15' 


The complete coefficient can thus be written as 


(4.135) 


CT = C'rP 
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1 


15 12 tt 4 6 12 J 4nb e 


dxj2 b - 
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n =14 


d n 
db r ’ 


6 6 arctan 


/ i m n T \ 

l ReflT J 


I b=§,y=l J 

' (4.136) 

This expression is in fact independent of b (see Eq. (4.134)). Therefore we may 
set simply 6 = 1. Summing up the (Taylor) series, we find after the substitution 

X = 47 TZ 3 



7r 

15 



/ 128 + 3tt 4 z 3 -8tt 2 (2 + 3z 2 ) 
1 67 t 3 z 


+z 2 arctan 


7r(l 

-2z + (z 2 - 



(4.137) 


From this expression we see that the complete HDL self energy is necessary for 
this coefficient (and not only the expansion for small go, which is sufficient for 
the fractional powers and the logarithmic terms). The remaining integral over 
the parameter z can probably not be done analytically. Numerically one finds 


c T = -0.00178674305 ... 


(4.138) 
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The computation of the longitudinal part proceeds in a similar way. With the 
substitution (4.79) the longitudinal gluon self energy can be written as 


= 5eff H 2 Hl 


(4.139) 


with some function Hl- In a similar way as in Eq. (4.134) we can make the 
estimate 

t? f (4140) 

. hr , \ x n — 3 


(4.140) 


In a similar way as above we obtain 


- _ -a) _ it _j_ r 

CL ° L 15 2tt 6 6 12 J nb 3 


, o / o , 

ax > — —— 6' 

4-' n! 56" 


x arctan 


Irnll j, 

Ilefr / 


I6=0,y=l / 

(4.141) 


4 = (^ 37r — 2 tt (12 — 7 T ) ^-17 + 67 e - 61og(vr) - 540 —j j . 

(4.142) 

We have to be careful when summing up the Taylor series in Eq. (4.141) for the 
following reason. The function RelR has a zero at 6 3 ~ 0.2653x, and therefore 
<t(6 ) := ar ctan (ImER, / Reft i) is not an analytic function of 6 . We have to replace 
cr( 6 ) with(f( 6 ) := — arctan(ReIlL/ImT[i / ) + f • cr( 6 ) and cr( 6 ) have the same Taylor 
expansion around 6 = 0, but only a is an analytic function of 6 . In this way we 
arrive at 


~ CL = ^-^l dZ 


ir(ir 2 + 12(—1 — z 2 + z 3 )) 
24i 


2 2 z l/^ + l 

+z arctan-1-log - 

7T 7T V 2 — 1 


~ 0.11902569216... (4.143) 


For the complete T 4 coefficient we also have to take into account the 0{T 4 ) 
contribution from Eq. (4.121). Thus we get finally 


C = C T + C L -- 

_ 2048 - 256?r 2 - 36 tt 4 + 3tt 6 
~~ 2160vr 2 


C'(4) 
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7 r 4 
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25920tt 2 + 2160tt 2 


-/ dz 
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1024 + 7 T 6 - 64tt 2 (2 + 3 z 2 ) - 12tt 4 (1 + z 2 - 3 z 3 ) 
247T 3 Z 
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+2z z arctan 


7 r 


(i-^) 


_2z + (z 2 -l)log(§±±) 
~ 0.0075766779858 ... 


+ z arctan 


2 z 1 fz +1 

-h - log -- 

7T 7T \ 2 — 1 


(4.144) 


4.5.8 How to compute higher order coefficients 

Using the methods described above one can compute the coefficients of higher 
terms in the expansion of C v with respect to b. This is straightforward for the 
coefficients of the fractional powers and the logarithmic terms, where one only 
has to include higher orders in the expansion of the HDL gluon self energy with 
respect to qo, see Eqs. (4.70), (4.71), (4.77) and (4.78). In this way one finds for 
the coefficients of b 14 , b 16 and 6 18 log b in the pressure 


C 14 = 


Cl6 — 




^ 18 ,log — 


128 2 2 / 3 (—70 + 97t 2 )T (M) C (^) 
243v/3vr 17 / 3 


~ 0.21345... , 


8 2 1 / 3 (197120 - 28800vr 2 + 729 tt 4 )T (f) C (f) 
3645 v/3tt 19 / 3 


(4.145) 


-0.75413 ... , (4.146) 

-94371840 + 15728640vr 2 - 606208vr 4 - 7920 tt 8 + 2160vr 10 - 135vr 12 

22680tt 2 


18.61898... . 


(4.147) 


We note in passing that the origin of the b 6 log b term is quite different from 
the origin of the b 6n log b terms with n > 1: the former arises from dynamical 
screening, while all the remaining logarithmic terms come from integrals of the 


dx 


type L . 

For the terms of the order b 12+6r , with r G No one has to sum up again IR en¬ 
hanced contributions, in a manner similar to the evaluation of the 6 12 -coefficient 
described above. In order to see how these IR enhanced terms come about, let 
us make again an estimate of the type (4.134), but now taking into account also 
the explicit q 2 and q$ from the inverse free propagator in Eq. (4.67). Using the 
fact that q/(g e sg) oc te 1 / 3 and qo/(g e tfl-t) oc b 3 y, we find in the transverse sector 
contributions of the type 


dx 


* 'E'min.T 


b 2 y 


x 


1/3 


^i/spy y)2 , 

ti 2p 3 


(4.148) 


When the denominator on the right hand side vanishes, the integrand on the left 
hand side is equal to fe 12 + 2 (p+ s )y 3 + 2 (P+ s ) / x, which gives logarithmic terms upon 
integration with respect to x (which are not IR enhanced). The important point 
is, however, that arbitrarily high powers of b in the integrand produce terms 
of the order ft 12 + 6r , where r is the total number of q 2 and q 2 t insertions in the 
integrand. 
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For the longitudinal part, the situation is completely analogous, as can be 
seen from the estimate 


' 'Emin.L 


dxx 2 (^pj x 2p {b 3 y) 2s 


foV2+G(p+s) y3+2(p+s) 

n — 2p — 3 


(4.149) 


As an illustration let us compute the coefficient of 6 18 . From Eq. (4.148) we 
see that we get contributions from (p. s) = (1,0) and (p. s ) = (0,1). Instead of 
simply setting the explicit q 2 and q 2 to zero as in Eq. (4.136), we have to consider 
now the object 
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Imn T (g 0 , g) 


g 2 + go + R eU T {q 0 ,q ) 


q=g o=o 


_ = (g 2 -go)Imnr(g 0 ,g) 

Imn r (g 0 , q ) 2 + ReIXr(go, g) 2 

Therefore we obtain in place of Eq. (4.136) 
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6=0,2/=! . 


where 
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(4.152) 


18 ’ T 127575vr 2 

-8(4158720 - 715680vr 2 + 29921 tt 4 ) + 675?r 6 


In a similar way as in Eq. (4.137), we can write Eq. (4.151) explicitly as 
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7T 2 (z 2 — l) 2 + ^2z + (z 2 — 1) log ^ - 

Numerically, one finds readily 

cis T = 7.97457343372231.... 


(4.153) 


(4.154) 
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Thus we obtain 

0^4 roo 2 

cis,l = cYs tL + — J dz - Y^(176 - 487r 2 + 37r 4 + 962; 2 

-12t rV + 48 z 4 ) +-—- 2 

7T 2 + (2® — log (f=t)) ‘ J 

~ 8.59505395976.... (4.157) 

To summarize, the temperature dependent part of the pressure at small tem¬ 
perature can be expanded in a series with coefficients of the orders g 4 /i 4 b <>+ " 
and g 4 /i 4 & 6 ( 1+r d log b, with n E No- In this and the previous subsections we have 
described in detail how to compute the coefficients of this series. On the other 
hand, the computation of terms which are suppressed by explicit powers of g e g 
would be much more involved. In particular one would have to include also 
diagrams beyond the large -Nf formula (4.65). 


4.5.9 Specific heat 

Neglecting terms which are explicitly suppressed with powers of g e g, the pres¬ 
sure can be written as 


P = Pfree + <7 e V/ (~) + 0{ff V) (4.158) 

V fleff g J 

with some function /. The specific heat (4.11) can be written as a sum of two 
terms, C v = C v g + C Vt 2 , where 


and 


C 


_ m d 2 P 

Cv d ~ T Qjg 2 ’ 


(4.159) 

f d 2 P \ 2 

/d 2 p 

(4.160) 

v ' 2 ~ T \dgdTj / 
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Using the expression (4.158) for the pressure, and setting T = b 3 g e g/j,, we find 


l-'r),l,/ree)/(fl l efth') 

= b 3 f"(b 3 ) 

(4.161) 

~ C V ,2Jree )/(fleffb ) 

= 0 (g 2 s ). 

(4.162) 


and 


This means that beyond the free contribution C v ,2 only contains contributions 
which are suppressed with explicit powers of g e g. Within our accuracy we may 
therefore neglect C V: 2 in the interaction part of the specific heat. 

Adding the transverse cut contribution [Eq.(4.74)], the longitudinal cut con¬ 
tribution [Eq. (4.81)], the part of the pole contribution [Eq. (4.121)] which is not 
exponentially suppressed, and the non-n/, contribution [Eq. (4.131)], we arrive 
at our final result the the specific heat at low temperature 5 . 


r - C° 

'-'V '-'v 


N n 


9> 2 T 

367r 2 




2 2 /3 r (I 

-40 V3 


) C (I) T 5/3 (5eff/x) 4/3 + 560 2V3r ( 3 H ( 3°) 


+ 


27x/37T n /3 

2048 - 256?r 2 - 36 tt 4 + 3t r 6 3 
180 tt 2 


log 


81\/37r 13 / 3 
9eSl^\ 


T 7 / 3 (9ezrf /3 


T J 


+ c 


+0{T u/3 /{g eS g) 2/3 ) + 0(<&/x 2 T log T), 


(4.163) 


where the constant c is given by 


2160t r 2 c 

2048 - 256?r 2 - 36 t r 4 + 3 tt 6 


7 

12 


~ -0.7300145937831 


(4.164) 


with c given in Eq. (4.144). We note that for sufficiently small temperatures there 
is a significant excess of the specific heat over its ideal-gas value, whereas ordi¬ 
nary perturbation theory [114, 115] would have resulted in a low-temperature 
limit of C v /C® = 1 — 2 a s /it. 


4.5.10 HDL resummation 

As we have seen in the previous subsections, the nonanalytic terms in the low- 
temperature expansion of the entropy density are determined by EIDL contri¬ 
butions to the gluon self energy. Terms beyond the HDL approximation are 
relevant for contributions from hard momenta q ~ g, yielding a term of order 
g 2 ffH 2 T 2 in the temperature-dependent part of the pressure, as shown in Eqs. 
(4.110) and (4.113). In order to retain all contributions that are nonanalytic in T 

5 The terms in the first two lines of Eq. (4.163) beyond the leading logarithm have been calcu¬ 
lated for the first time by A. Ipp [77]. 
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Figure 4.3: The first few orders of the low-temperature series for the entropy density in 
comparison with the full HDL-resummed result (from [113]). 
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The sum of the transverse and longitudinal contribution decays like q 3 for large 
q, therefore we may send the upper integration limit of the (/-integration to infin¬ 
ity. This just amounts to dropping terms that are suppressed by explicit powers 
of ste¬ 
in Eq. (4.165) we have neglected terms like g^gb 6 = (T ///) 2 in P/(g* s y 4 ), 
therefore we have to assume T <C p. This is however a weaker condition for the 
temperature than the condition T <C g e sp which we had in the previous sub¬ 
sections. The weaker condition is sufficient now, because we do not perform an 
expansion in b in this subsection. Adding the soft contribution from Eq. (4.165), 
the hard contribution from Eqs. (4.110) and (4.113), and the non-n& contribution 
from Eq. (4.131), we arrive at the following result for the entropy density [113], 
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where 5° is the ideal-gas entropy density. The right hand side of Eq. (4.166) is es¬ 
sentially given by one universal function of the dimensionless variable T / (g e {ip), 
which we define through [113] 

N 8 f 2 T ( S - S °) =■■ © (—) + Otfes)- (4-167) 

N g gi s ^T \9esgJ 

and which we have normalized such that the ordinary perturbative two-loop 
result [114, 115] for the low-temperature entropy density corresponds to 6 = 
— 1. The function © can be evaluated numerically. Fig 4.3 shows a comparison 
of the first few orders of the low-temperature series with the full HDL result 
taken from [113]. Further numerical results, in particular a comparison with 
non-perturbative large-Ay results, can be found in [100,113]. 


4.6 Some remarks on neutrino emission from ungapped 
quark matter 


A young neutron star loses energy mainly via neutrino emission from the bulk 
[33,57]. One generally distinguishes fast processes, which lead to neutrino emis- 
sivities proportional to T 6 , and slow processes which lead to neutrino emissiv- 
ities proportional to T 8 , see e.g. [33] and references therein. The dominant pro¬ 
cess for neutrino emission from ungapped quark matter is the quark analogon 
of the direct Urea process [57,116,117], 


u + e 


u + e +u e , 
d -\- u e . 


(4.168) 

(4.169) 


The cooling behavior is governed by the cooling equation 

dT 

C V (T)— = —e(T), 
where e(T) is the neutrino emissivity, which is given by [57,116, 117] 

d 3 p d 1 f d 3 p u 1 f d 3 p e 1 f d 3 p u 1 
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xn f (E d - p d )[l - n f (E u - p u )][ 1 - n f (E e - g e )\ 

+ \Mec\ 2 (2lT) 4 5 4 (p u +p e -p d ~ p u ) 


(4.170) 


■Ev 


xjif(yE u g u )nf(E e /i e )[l nf(E d Pd)\ 


(4.171) 


where n t are the fermion helicities, E t are the energies of the fermionic quasipar¬ 
ticles as determined from the poles of the propagators (neglecting the imaginary 
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parts of the self energies), Mp is the matrix element for the reaction (4.168), and 
M ec is the matrix element for the reaction (4.169). 

In [57] it is shown that the emissivity through leading logarithmic order can 
be written as 


e(T) = 


2 g 2 Gp cos- 9 C 2 6 


7T 


/ oo rc 

dx d / 

-oo J — c 


dxn, 


dx u x^ 


xnf(Tx d )nf(-Tx u )rif(T(x u -x d + x u )) v l (E d ) v \E U ) 




,(4.172) 


where Gf is the Fermi coupling, 6 C is the Cabibbo angle, and \x q ~ ~ ft,/. 

The inverse group velocity v~ x , as defined in Eq. (3.103), contains non-Fermi- 
liquid contributions, which we have discussed in chapter 3. As shown in [57] 
the non-Fermi-liquid corrections to the specific heat and the emissivity lead to a 
(modest) reduction of the temperature at late times. 

If one wants to go beyond the leading logarithmic accuracy of Ref. [57], at 
least the following ingredients will be required. First, there will be corrections 
to the emissivity from gluonic corrections to the weak interaction vertex [57], It 
might also be necessary to generalize the formula (4.171) for the emissivity in 
a way that takes into account the finite lifetime of the fermionic quasiparticles, 
since the imaginary part of the quark self energy is of the same order as the term 
which fixes the scale under the logarithm in the real part. 

As a first step towards a more complete computation of the emissivity, let 
us take the scale under the leading logarithm of the group velocity as given by 
Eqs. (3.67) and (3.68), and simply insert the resulting expression for the group 
velocity into Eq. (4.172). Using the integral 
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(4.173) 


(4.174) 


While the constant under the logarithm will certainly be modified in a more 
complete computation as discussed above, it is already interesting to note that 
this constant is about an order of magnitude larger than the corresponding con¬ 
stant in the specific heat, which we find to be approximately equal to 0.282 from 
Eq. (4.163). 


66 



Chapter 5 

Color superconductivity 


5.1 Basics of superconductivity 

Since its experimental discovery in 1911, superconductivity has been extensively 
studied both in experimental and theoretical condensed matter physics. The first 
microscopic theory of superconductivity was developed by Bardeen, Cooper 
and Schrieffer (BCS) in 1957 [118], see e.g. [106, 119, 120] for pedagogical re¬ 
views. The key observation is that the electron-phonon interaction induces an 
effective electron-electron interaction which is attractive 1 . This leads to the for¬ 
mation of Cooper pairs. The ground state, which is different from the ground 
state of a normal Fermi liquid, is then given by 

i^)= n ( Wk + wkc kT c -kj.) i°)> t 5 - 1 ) 

k 

where pL and cj^ are creation operators for electrons with spin up and spin 
down, respectively. The coefficients Uk and t>k are determined by minimization 
of the free energy. The state (5.1) can be characterized by a non-vanishing "di¬ 
electron" condensate, 

W) + 0, (5.2) 

which leads to a gap in the fermionic quasiparticle spectrum. A condensate 
like (5.2) breaks spontaneously the 17(1) invariance of electrodynamics. There¬ 
fore the photon will acquire a mass via the Higgs mechanism. This leads to the 
Meissner effect, namely the screening of magnetic fields inside a superconduc¬ 
tor. 

1 Using Wilsonian renormalization group methods [121], it can be shown that the four fermion 
operator corresponding to the scattering of two electrons with opposite momenta is the only 
marginal operator in the vicinity of the Fermi surface, all other fermionic interactions being ir¬ 
relevant [122]. 
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5.2 Basics of color superconductivity 

Let us begin with an ensemble of free quarks. At finite density the ground state 
of this system is simply given by a Fermi sphere. Now let us turn on inter¬ 
actions. At asymptotic densities at first sight a purely perturbative treatment 
might seem appropriate because the interactions are "weak" due to asymptotic 
freedom. However, it is important to note that single gluon exchange is attrac¬ 
tive in the color anti triplet channel 2 . This has dramatic consequences, since an 
arbitrarily weak attractive interaction leads to the Cooper instability of the Fermi 
surface, and the ground state will contain diquark condensates, 

p)y^(-p)} ± 0, (5.4) 

where i, j are fundamental color indices, /, g are flavor indices, a, (3 are spinor 
indices, and 'ib c denotes the charge conjugated spinor, L c = C'h T . By virtue of 
its similarity to superconductivity in ordinary condensed matter physics, this 
phenomenon was termed color superconductivity. 

In the usual notation for the irreducible SU(3) representations one has 

3 <g> 3 = 6 © 3. (5.5) 

Therefore one cannot construct a color singlet from the expectation value (5.4). 
This means that global color symmetry is spontaneously broken by the diquark 
condensate. 

Quarks have color, flavor and spin as internal degrees of freedom. Therefore 
quite many different color superconducting phases could be possible, depend¬ 
ing on the color, flavor and spin structure of the expectation value (5.4). At very 
high densities the phase with the lowest free energy is the so-called color fla¬ 
vor locked (CFL) phase [123,124,125,126], where up, down and strange quarks 
contribute to the pairing on an equal footing, 

$CFL OC £ ijA £ f9A = 6 if S j9 - 5 i9 6 jf . (5.6) 

The symmetry breaking pattern of the CFL phase is (neglecting electromag¬ 
netism and quark masses for the moment) 

SU (3) c x SU(3) l x SU(3) r x U(l) B -> SU(3) C+L+R x Z 2 . (5.7) 

Here SU( 3) c corresponds to global color symmetry, SU(3) Rj l correspond to 
right and left handed flavor symmetry, and U{l) R corresponds to baryon num¬ 
ber conservation. These symmetries are broken down to the diagonal subgroup 

2 This can be seen as follows. The tree level scattering amplitude of two quarks is equal to the 
tree level scattering amplitude of two electrons in QED, times the QCD factor 

TikTji = — — (5ikSji — SuSjk) + g {fiikSji + Su5jk). (5.3) 

Since the electron-electron interaction is repulsive, the negative sign of the first term in Eq. (5.3) 
implies that the quark-quark interaction is attractive in the antitriplet channel. 
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SU(3) c+ l+r times a 77 group that corresponds to vj —ip. In the CFL phase 
all eight gluons acquire a mass through the Higgs mechanism [127, 128, 129]. 
The electromagnetic U{ 1) symmetry is also broken, but there remains an un¬ 
broken t7(l) [123, 130], whose generator is given by Qcfl = Q + -^Tg, with 
Q = diag(— |) in flavor space, and T 8 = -^=diag(l, 1, —2) in color space. 

As in every field theoretical model with spontaneous symmetry breaking, 
Nambu-Goldstone bosons appear in the CFL phase [127,131,132,133,134,135]. 
Because of the symmetry breaking pattern (5.7) there are 17 broken genera¬ 
tors, which would lead to 17 Nambu-Goldstone bosons. Eight of these, corre¬ 
sponding to the breaking of SU (3) c , are "eaten" by the gluons. The "surviving" 
Nambu-Goldstone bosons are an octet corresponding to flavor symmetry break¬ 
ing and a singlet corresponding to baryon number breaking. While this singlet 
boson is strictly massless, the octet bosons become massive if finite quark masses 
are taken into account. Furthermore there is a light singlet from 17(1 )a breaking. 

At lower densities, the strange quark will eventually decouple, which should 
lead to a two-flavor color superconducting (2SC) phase [38, 39, 40], where only 
up and down quarks participate in the pairing, 

4>2sc oc £*J 3 e /s ; y (5.8) 

In the 2SC phase the global color symmetry group SU (3) c is broken down to 
SU (2) c , but flavor symmetry remains unbroken. Therefore there are five would- 
be Nambu-Goldstone bosons [136, 137, 138], which are "eaten" by the gluons. 
The gluons with adjoint colors 1-3 are massless, while the gluons with adjoint 
colors 4-8 acquire a mass through the Higgs mechanism [139,140]. As in the CFL 
phase the electromagnetic 1/(1) symmetry is broken, but again an unbroken U( 1) 
survives [136,141], whose generator is now given by Q 2 SC = yyMl — 2\f3Ts). 

In Ref. [142] it was argued that the 2SC phase will be disfavored if one takes 
into account the constraints from color and electric charge neutrality at finite 
strange quark mass m s . To date a completely reliable description of the phase 
structure of QCD at densities below the CFL regime is still lacking. In the fol¬ 
lowing we will list some of the phases that have been proposed in the literature. 

Computing the meson masses from the low energy effective Lagrangian in 
the CFL phase, one finds that the masses of the K + and K° become imaginary 
if m s exceeds a certain value. This indicates the formation of a kaon condensate 
[44,143,144,145]. 

The charge neutrality condition and the finite strange quark mass lead to 
a mismatch of the Fermi momenta of the different quark flavors. This might 
give rise to so-called gapless superconducting phases (g2SC [47, 48, 146, 147] 
and gCFL [148, 149, 150, 151]), where the fermionic quasiparticle dispersion 
laws are modified, some of them corresponding to gapless excitations. In Refs. 
[152, 153, 154] however, it was found that the Meissner masses are imaginary 
in these phases. This would indicate that the gapless phases are unstable, if no 
mechanism can be found which makes the Meissner masses real quantities. 
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So far we have assumed only spin zero diquark condensates. The gaps re¬ 
sulting from spin one condensates are much smaller than in the spin zero case. 
Nevertheless spin one condensates might play a certain role at intermediate 
densities, since they may be formed from quarks of the same flavor, see Refs. 
[38,155,156,157,158,159]. 

The mismatch of the Fermi momenta might also induce a so-called LOFF 
(Larkin-Ovchinnikov-Fulde-Ferrell) pairing [56, 160, 161, 162, 163]. The LOFF 
state is characterized by Cooper pairs with non-zero total momentum, and there¬ 
fore translational and rotational symmetry are spontaneously broken. This leads 
to crystalline structures that might have relevant consequences for astrophysics 

[41]. 


5.3 The color superconductivity gap equation 


The gap equation for a color superconducting phase can be derived from the 
Schwinger-Dyson equation for the quark two-point function, allowing for a non¬ 
vanishing expectation value (5.4). At this point it is convenient to introduce a 
Nambu-Gor'kov basis for the quark spinors, which is given by T = (L, 'il' c )' r , 
T = (L, L c ). The inverse quark propagator in this basis is given by [38,164,165] 

Q + M7o + 

<&+ $- 



7o + S 


(5.9) 


where $ ± are the gap functions, related by 4> _ (Q) = 7 o[$ + (Q)]^ 70 / and E(Q) 
is the quark self energy, with E(Q) = G[E(—Q)] r G _1 . Flavor and fundamental 
color indices are suppressed in Eq. (5.9). 

In order to solve the Schwinger-Dyson equation a suitable approximation 
scheme has to be chosen. To date only the exponent and the leading contribu¬ 
tion to the prefactor of the gap [see Eq. (5.13) below] have been determined 
consistently, where the following set of approximations turned out to suffice: 

• replace the full quark-gluon vertex with its tree level counterpart (at least 
in Coulomb gauge) [166], 

• replace the full gluon propagator with its HDL approximation [167], 

• replace the quark self energy with the leading logarithmic result of the 
normal phase (see Eq. (3.67)) [164], 

A rather lengthy calculation [45,156,164] then yields the following gap equation 
for the 2SC phase (neglecting the gauge dependent part) 


(j) + (e k ,k) 


18vr 2 


f d(q — n)Z(eq )—— ^ q ’ ^ tanh 
Jo e q 




(5.10) 
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where ek 
and 


y 7 (k — g) 2 + \4> + \ 2 , 5 ~ 0(gn ) is a cutoff, b = 2567r 4 (2 /{Nfg 2 )) 5 / 2 , 


Z(E — g) 


dRe£+ 

9(^-/x) 


(5.11) 


is usually called "wave function renormalization factor" in this context (in the 
previous chapters we referred to this quantity as group velocity). With the defi¬ 
nition 


S- log ( ^ f 

3\/27r \k — g + £k / 


(5.12) 


the solution of the gap equation at zero temperature can be written as [86, 156, 
164,168,169,170,171] 


<t> + 


2bb' 0 gexp 



(5.13) 


where b' 0 = exp[— (n 2 + 4)/8]. F(x) is quite a complicated function, which is 
given explicitly in [164]. At the Fermi surface F(x) is equal to 1 + 0(g 2 ), and 
away from the Fermi surface it decreases rapidly. 

The parametric dependence of the gap on the coupling constant as shown 
in Eq. (5.13) was first derived by Son [168] using renormalization group tech¬ 
niques. It is different from the BCS result <j) + ~ exp(—1 /g 2 ), which is a conse¬ 
quence of the fact that quasistatic chromomagnetic interactions are only dynam¬ 
ically screened. 

The result for the gap can also be derived within the framework of high 
density effective theory (HDET), in which the only relevant fermionic degrees 
of freedom are those in the vicinity of the Fermi surface [43,172,173,174, 175]. 

In the literature also the possible gauge dependence of the prefactor of the 
gap, which would indicate an inconsistent approximation scheme, has been dis¬ 
cussed. In [176] the gap equation was solved numerically in a general covariant 
gauge, with the result that large gauge dependences occur for g > g c ~ 0.8. In 
Ref. [177], however, it was argued that in a covariant gauge one has to include 
the quark gluon vertex correction in order to obtain a gauge independent result. 
In [178] it was shown that the prefactor of the gap is gauge independent in a 
general Coulomb gauge, if the gap is evaluated on the quasiparticle mass shell. 
In the next section we will give a general proof that the fermionic quasiparticle 
dispersion laws in a color superconductor are gauge independent [179]. 


5.4 Gauge independence of fermionic dispersion laws 


The inverse quark propagator in the Nambu-Gor'kov basis, as shown in Eq. 
(5.9), is the momentum space version of the second derivative of the effective 
action. 


<5 2 r 




i/,=y=A“=0,Ag=Ag’ 


(5.14) 
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where T = (ip, ip c ) T , T = (ip, ij; c ), and Aq is the expectation value of Aq, which 
is in general non-vanishing (see Sec. 5.6). 

It should be noted that the doubling of fermionic fields in terms of T and ^ 
is just a notational convenience here; the effective action itself should be viewed 
as depending only on either (ip, ip) or the set 'E = (ip, ip c ) T . 

The gauge dependence identity for the effective action [Eq. (2.12)] can be 
written as 


ir = I d ‘ x {sW) sx w {x) ~ sx (*^ x) sW) + sJ^) Sx w (x) ) 

= f* x (s^x) SX ^ (x) + + ' <5 ' 15) 

Eq. (5.15) can be cast in a more compact form using the DeWitt notation, but 
now only for the fermions, 

«r = r,W + J (5.16) 


where i = (ip(x), ip c (x)) T , i = (ip(x),ip c (x)), and the comma denotes functional 
derivation. Taking the second derivative of (5.16), setting ip = ip = Af = 0, 
Aq = .4[j, and using the fact that 


ST 

Mg 


i/,=i/,= J 4g=0,Ag=Ag 


= 0 , 


(5.17) 


we obtain a gauge dependence identity for the inverse propagator (5.14), 


STij = -r k 3 sx* 


+ T ki 6Xj + 
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d A x 


ST 
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SA a0 (x 


- 8 X$)(x)- 


(5.18) 


Up to this point, our functional relations are completely general and apply 
also to the case of inhomogeneous color superconducting phases. We do not 
attempt to cover the complications this case may add to the question of gauge 
independence, but continue by assuming translational invariance. This allows 
us to introduce 5A a0 := —6X( ,( P (x = 0) and to write (5.18) as 

dT - 

5T fj = -rp-SX^ - 5X k f ik - gj^ 6 A a0 . (5.19) 

Furthermore, we can transform Eq. (5.19) into momentum space, 

5T fj (Q) + SA a0 ^'0>- = —T k j(Q)5X k i (Q) - 8XK(Q)T ik (Q), (5.20) 

where the indices i and i from now on comprise only color, flavor, Dirac and 
Nambu-Gor'kov indices. Using Eq. (C.l) we obtain 

6 det(Tjj) + SA aQ det(Tjj) <5 tot det(T iJ ) = - det(T^)[^ + ^]. (5.21) 
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The left hand side of this identity is the total variation [62, 63] of the determinant 
of the inverse quark propagator, with the first term corresponding to the explicit 
variation of the gauge fixing function, and the second term coming from the 
gauge dependence of Aq. 

Since the determinant is equal to the product of the eigenvalues, Eq. (5.21) 
implies that the location of the singularities of the quark propagator is gauge 
independent, provided that the singularities of 5X k k do not coincide with those 
of the quark propagator. The singularity structure of SX may be discussed in a 
similar way as in the case of the Higgs model, see Sec. 2.3.3. The quantity SX 
is 1PI up to a full ghost propagator, and up to gluon tadpole insertions (see Ref. 
[60] for the explicit diagrammatic structure) 3 . As in [60, 64] one may argue that 
the singularities of the ghost propagator are not correlated to the singularities of 
the quark propagator. Gauge independence of the zeros of the inverse fermion 
propagator then follows provided that also the 1PI parts of SX have no singular¬ 
ities coinciding with the singularities of the fermion propagator. An important 
caveat in fact comes from massless poles in the unphysical degrees of freedom 
of the gauge boson propagator, which are typical in covariant gauges and which 
can give rise to spurious mass shell singularities as encountered in the case of 
hot QCD in [73]. But, as was pointed out in [74], these apparent gauge depen¬ 
dences are avoided if the quasiparticle mass-shell is approached with a general 
infrared cut-off such as finite volume, and this cut-off lifted only in the end, i.e., 
after the mass-shell limit has been taken. 

The determinant which appears in (5.21) is taken with respect to color, flavor, 
Dirac and Nambu-Gor'kov indices. The determinant in Nambu-Gor'kov space 
may be evaluated explicitly using Eq. (C.2). From Eq. (5.9) we obtain in this 
way 

det(T-) = det[($-p7o + S)($ + p7o + I!) 

-(<$ ~ P7o + £)$"(# - P7o + S)- 1 4>+], (5.22) 

The inverse of the matrix of which the determinant is taken here appears, of 
course, in the ordinary quark propagator, which is obtained by inverting (5.9), 

G + (Q) = [{Q — P7o + E)($ + pyo + E) 

~(Q - ^70 + £)<&"($ - Mo + sr^+r 1 ^ - ^7o + S). (5.23) 


At leading order, when the quark self energy E can be neglected, (5.22) can 
be approximated by 

det(T-) ~ det [(Q - p7o)($ + P7o) ~ {Q ~ Mo)$~ - /x7o) _1 $ + ] • (5.24) 


This expression can be rewritten as 


det(Tjj) ~ det (qjj - (g - p) 2 - (0+ + )V+ + ) 


V 


r+ 


3 It should be noted that SX k k will not contain any gluon tadpoles if its perturbative expansion 
is constructed using the shifted field A' = A — A. 
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(5.25) 


+ (?o - (? - ^) 2 - (^) t </> z + _) 

+ (?0 — (^ + ^) 2 ~~ (^r-)Vr-) ^V- 

+ (<7o _ (9 + m) 2 “ (<^+)Vz+) ^z + 

where the projection operators introduced in [180] are given by 

V% ± =Vr,iT ± ( q)A ± (q), (5.26) 

where V r ,i = 5(1 ± 75) is the chirality projector, P±(q) = |(1 ± 75707*9*) is the 
helicity projector, and A =b (q) is the energy projector given in Eq. (3.42). (Note 
that V+_ = V+ + = V~ + = V^_ = 0 in the massless limit which we are considering 
here [180].) If we assume for simplicity that (jrcj) is diagonal with respect to color 
and flavor indices, the determinant in Eq. (5.25) can be readily evaluated, 

det(Ejj) ^ n^-^H^I 2 ) (,o 2 -(,-.) 2 -|^ (/) | 2 ) 

/ 

X (<7o - (<? + V? - |0r- f) | 2 ) (90 - (? + m) 2 - l^+ (/) | 2 ) , (5-27) 

giving the well known branches of the quasiparticle excitation spectrum for each 
value of the color-flavor index / [180]. 

At leading order, when the quark self energy can be neglected, Eqs. (5.21) 
and (5.22) imply that the gap function is gauge independent on the quasiparticle 
mass shell. It should be noted, however, that at higher orders only dispersion 
relations obtained from (5.22), which also include the quark self energy E, can 
be expected to be gauge independent. 

We note furthermore that for 4> ± —► 0 the determinant (5.22) becomes 

det [(Q + P70 + £)($ - ^70 + E)] = det(^(Q) _1 ) x det(^(-Q) _1 ), (5.28) 

where Q {) ( Q ) = [(5 + M7o + E] ~ 1 . Therefore, in this case the gauge independence 
proof, of course, boils down to the gauge independence proof for the singulari¬ 
ties of the quark propagator without diquark condensates. 


5.5 Computation of 8X 

The aim of this subsection is the evaluation of SX k k for the 2SC phase in covari¬ 
ant gauge. This is interesting because the existing calculations in the literature 
give gauge dependent results for the gap in this gauge 4 . 

4 In [181] the gap equation was considered in a non-local gauge, with the result that the "best" 
value of the gauge parameter in covariant gauge should be £ ~ |, but then the result for the gap 
does not agree with the well-established result in Coulomb gauge [178]. 
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Figure 5.1: Feynman diagram for SXn^ (see Eq. (5.29)). 


We will focus on possible singularities of 8X k k , since they might invalidate 


the above proof of gauge independence. At one-loop level 8X k k contains terms 
like 


5X W (K) := g 2 


+T 


d 4 Q 
i(2ir) 4 
1 1 


Tv D,c,f 


T n l,^P^D^(P)G + (Q)YT a 


a P 2 s/2f 


P' i D ia ,(P)Q-{Q)'fT'Z 


P 2 V2f 


(5.29) 


where P = K — Q, l) jW is the gluon propagator, G~ is the quark propagator (see 
below), and we have used 


D l 

s->a 

yp 

5F a oc F a 


r 5 ab , 


-igT a lf>, j P T a V’c 

1 

P 2 ' 

\/2f ^ 


(5.30) 

(5.31) 

(5.32) 


Eq. (5.29) corresponds to the diagram in Fig. 5.1, where the vertex with a circle 
denotes 5F a and the small triangle denotes D l a . 

In the 2SC phase the gap matrix is given by [39,156] 

l* + )l 9 = e^e ij3 (f + (Vf + -V+) + - Vf + )), (5.33) 

where /, g are flavor indices and i, j are fundamental color indices. The Vs are 
the projection operators defined in Eq. (5.26). We have assumed for simplicity 
that the right-handed and the left-handed gap functions are equal up to a sign 
[156]. For the quark propagator G ± {Q ) one finds then [139] 


with [139] 


[OH 9 = $ f9 [(<% " ^p)G ± + 5 i3 8 j3 G±] , 


G ± (Q) 

Gq(Q) 


E 

e=d= 

E 

e=d= 


go T (h - 

g 0 2 -(h-eg ) 2 

go T (h ~ eg) 

go (m eg) 2 


£ q^2 A±e(q)70 ’ 

A±e (q)7o, 


(5.34) 


(5.35) 

(5.36) 
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where are the energy projectors given in Eq. (3.42). We shall consider only 
the part of <5A' (1 ) which contains the propagator G ± , because it is this part which 
depends on the condensate. We will call this part JApp. In the HDL approxi¬ 
mation the gluon propagator obeys the following Ward identity. 


P^D^P) 


’ p '2 


(5.37) 


Thus we obtain after taking the color trace, neglecting the antiparticle propaga¬ 
tor and the imaginary part of the gap function [156], and dropping an irrelevant 
prefactor 


SX (n) (K) = 2 g 2 



d 4 Q 1 _(Zo_ 

i(27r) 4 (P 2 ) 2 q ) 2 


^pTVoIA+707 v\P v 


WI 


\[i 


d 4 Q 1 _ qo(po ~ P • q) = 

i(27r) 4 (p 2 - p 2 ) 2 q 2 - (p - q) 2 - (</>+) 2 
d 4 Q — iq 4 (i(k 4 — q 4 ) — kt + q) 

i(2n) 4 [(k 4 - q 4 ) 2 + ( k 2 + q 2 - 2kqt)] 2 (qj + (p - q) 2 + {(j) + ) 2 ) ’ 

(5.38) 


where t = cos 6 and q 4 = —iqo- The integral is dominated by 0 ~ 0 , therefore 
we may set t = 1 in the numerator. We set k = // and define ( = q — //. After 
carrying out the (^-integration (assuming that (l> + is almost independent of q 4 ) 
we perform the analytical continuation k 4 —> iko and take the on-shell limit, 
A'o —> . Then the ^-integral (with UV cutoff ±//) gives terms proportional to 
log(A;Q — (4> + ) 2 ). Thus there may be indeed a mass shell singularity in 5X, at least 
within our approximations. It is however only a logarithmic singularity, which 
is harmless because on the right hand side of the gauge dependence identity we 
have expressions like 

{kl - (0 + ) 2 ) (log(fcg - (</> + ) 2 ) + finite) , (5.39) 

which vanish for k$ —> 7 . We conclude that no IR cutoff should be necessary to 
assure gauge independence of the dispersion laws (at zero temperature). 

An explicit calculation shows indeed that a simple IR cutoff for the gauge 
dependent part of the gluon propagator in the gap equation is not sufficient to 
remove the gauge dependence of the gap in covariant gauge. In Ref. [177] it has 
been argued that this gauge dependence will only disappear after the inclusion 
of vertex corrections in the gap equation. 


5.6 Gluon tadpoles 

In a color superconductor global color symmetry is dynamically broken by the 
diquark condensate. Therefore there is no symmetry which forbids the existence 
of gluon tadpoles, cf. the discussion in Sec. 2.2. This section is devoted to a 
detailed analysis of the tadpole diagrams. 
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Figure 5.2: Leading order tadpole diagram. 


5.6.1 Leading order tadpole diagram (2SC) 

Let us consider the one-loop tadpole diagram as shown in Fig. 5.2, 

T a = -f J ^4TVD,c,/,iVG[fSS'(Q)]. (5.40) 

Here S(Q) is the quark propagator in the Nambu-Gor'kov basis, which is given 
by Eq. (5.9). At leading order we may neglect the quark self energy £. The 
quark-gluon vertex FJj is given by 



r a 
1 o 



(5.41) 


with F") = 7 oT" and fg = — 7o(T ") 7 . The trace in (5.40) has to be taken with 
respect to Dirac, color, flavor and Nambu-Gor'kov indices. First we evaluate 
the trace over Nambu-Gor'kov space which gives 

ja =-\S i §^ Ti ' D ’ cj[ra ° G+{Q)+(5 - 42) 

where is given in Eq. (5.34) for the 2SC phase. Evaluating the trace over 
flavor and color space we get 

T a = g(T a ) 33 j ^ i Tr D [ 70 (G + - G~ - G+ + )]. (5-43) 


Assuming that 4> ~ 0 [156] and that b 4 has negligible four-momentum depen¬ 
dence in the vicinity of the quasiparticle pole we obtain 


T a ~ -4 g(T a 


133 


p-q 


p-q 


-^( T a ) 33 [ dqq 2 {p-q) 
7T 2 Jo 


d A Q 

*( 2? 0 4 V?o “ (A* “ d) 2 “ l^ + l 2 do - (F “ d) 2 

1 1 


Viv-d) 2 + \0 + \ 2 \d-d\ 


(5.44) 


In order to obtain an order of magnitude estimate, we make the approximation 
4> + (q) — 4 >q&( 2p — q) with = const. [139]. Then the ^-integration can be 
readily performed, with the result 


T a ~ -%{T a ) 33 

7T Z 


t* (^) 2 log (0) + 0((gp^) 2 ), (T a ) 33 = - 6 a8 /V3. 


(5.45) 












This result is in fact of order //((/>[j") 2 , because log(^/(2/x)) is of order 1/g [see 
Eq. (5.13)]. 

5.6.2 Gluon tadpoles for CFL with m s ^ 0 

We would like to evaluate the gluon tadpole in the CFL phase with a non¬ 
vanishing strange quark mass m s . The inverse Nambu-Gor'kov propagator 
takes the same form as in Eq. (5.9) (with £ —> 0 at leading order), where the 
free inverse propagators are now given by 

([fitf]" 1 )# = M + VloW s - S r3 S s3 ) + {0 + /i 7 o - m s )6 r3 S s3 ]8 ij , (5.46) 

([^] _1 )p = m ~ WoW s - S r3 5 s3 ) + (<%- P7o - m a )<r 3 $ s %, (5.47) 

Again the lower indices are color indices and the upper ones are flavor indices. 
At finite rn s the Dirac structure of the gap is more complicated than in the mass¬ 
less case [180], but for m s <C // we may assume that the additional terms are 
suppressed at least with m s /n [182], 

[* + ]£ ~ ^.+A+ 75 (SIS] - STS?) + O . (5.48) 

As usual we have 4>~ = 7 of 4> 1 ]^ 7 o- The inverse full quark propagator can be 
written as 


[g+]^ = [ g+)-w^-g-^ 


( bl 

h 


b 2 


b 2 



\ b % 63 


h \ 



h 


. (5.49) 


67 / 


The nine rows and columns of this matrix correspond to color and flavor indices, 
namely (color, flavor)=(l,l),(l,2),(l,3),(2,l),(2,2),(2,3),(3,l),(3,2), (3,3). The 6* are 
matrices in Dirac space, given explicitly by 

bi = ^ + P7o + 7o^ t 7o[^-P7o]^ 1 0 + 7o0 t 7o[^-P7o-«A]“V, 
b 2 = 7o<^7o[$ - Wo ~ m s ]~ l 4 >, 

6 3 = 7 o^ 7 o [0 ~ M 7 o] _ V, 

64 = 0 + P70 + 70^70 [0 ~ P 7 opV, 

65 = ^ + P7o-"is + 7o0 t 7o[<^-P7o]” 1 0 ) 
be = $ + /Do + 7o<p7o[$-/Do -m s ] _ V, 

h = 0 + P7o ~m s + 270 ^ 70 [#-/x7o] _ V, (5.50) 
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with 4 > = oi + A+75. The full quark propagator is obtained by inverting the matrix 
in (5.49), 


G + = 


( <J \ 0,2 03 \ 

O 5 

Ofi 

O 5 

02 a\ (13 

Ofi 

a 7 

07 

(i 4 (J 4 (lg J 


(5.51) 


with 


a\ = {bi - b 2 )~ 1 + a 2 , 

o 2 = (-263 + b 7 bz 1 (bi + b 2 )y 1 (b 3 - b 7 bz 1 b 2 )(b 1 - 6 2 ) _ | 

03 = (2f»3-M3 1 (^i + ^))” 1 , 

a 4 = -^(M&i -b 2 y 1 + (61 + 62)^2), 

05 = b q 1 , 

a 6 = 65 \ 

a 7 = ftg 1 , 

a 8 = -fry 1 ( 6 i + 6 2 )a 3 . (5.52) 


Similar expressions can be obtained for G := [[G Q ] 1 — 4> G ( \ 4> ] 1 . 

The gluon tadpole is given by 

T*:--§/([6 + ]"(T“), i - [6-]r/(n«) «”• (5.53) 


Taking the trace over color and flavor space we find that the only non-vanishing 
tadpole is T 8 , 


T 8 


9 

2\/3 


d 4 Q 

z(27t) 4 


Tr D 7o(ai + a 5 + a 6 


2a 7 — as) — (similar terms from G ). 


(5.54) 

We perform a Taylor expansion in m s , keeping only the leading term which is of 
order m 2 s . With the approximation cl) ~ 0 we get 


T a ~ 5' 


a8 


2v/3vr 3 J- 00 


dq4 / dg 


_ Q (0 + ) 2 (g| ~ % ~ d) 2 ~ 8(0 + ) 2 ) _ 

(<?4 2 + (9 - 9? + (0 + ) 2 ) 2 (d 4 2 + (q - d) 2 + 4(0+) 2 ) 


(5.55) 


This integral seems to be proportional to fp + ) 2 . It turns out, however, that the 
factor (d + ) 2 cancels against a factor l/(d + ) 2 from the result of the integration. 
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This is related to the fact that for (j) + —► 0 the integral would contain an infrared 
singularity, which is regulated by the gap [183]. 

We make the approximation p + — <pg 0(2// — q). After the integration we 
take the limit —> 0, since we assume m s S> (po- This leads to [184] 

T 8 ~ -nm 2 a 9^2 ( 21 “ 81o § 2 )- ( 5 - 56 ) 

5.6.3 Color neutrality 

To address the question of color neutrality [142,146, 185, 186, 187, 188] we con¬ 
sider the partition function 

exp(—fi/T) = J T>(.p exp(—5[<p]), (5.57) 

where tp denotes the set of all fields, and ,S'[y] is the QCD action (including gauge 
fixing terms and ghosts). Following the argument given in [189] it is easy to see 
that the system described by this partition function is color neutral, at least if 
one chooses a gauge fixing which does not involve Ag, for instance Coulomb 
gauge: The fields Ag appear in the action as Lagrange multipliers for the Gauss 
law constraint [190]. Therefore in the path integral the integration over the zero- 
momentum modes Ag 0 produces delta functions, S(N a ), where N a are the 
color charges. This means that only color neutral field configurations contribute 
to the partition function, q.e.d. 

The gluon tadpole in the 2SC phase [Eq. (5.45)] induces a non-vanishing 
expectation value for the gluon field Ag, 

A° ~ m~ D 2 ^ + f ~ (</> + ) 2 /(<?V), (2SC) (5.58) 

where is the corresponding leading-order Debye mass [139]. This expecta¬ 
tion value acts as an effective chemical potential for the color number 8 , 

li 8 = gA° 8 ~(cf> + ) 2 /(gii). ( 2 SC) (5.59) 

It may be noted that the chemical potential which has been found by requir¬ 
ing color neutrality in an NJL model is also proportional to </> 2 [187], In a similar 
way one finds in the CFL phase from Eq. (5.56) 

p 8 ~ m 2 Jn, (CFL) (5.60) 

which agrees again qualitatively with the result of the NJL model calculation 
[187], One should emphasize, however, that whereas in NJL models color neu¬ 
trality has to be imposed as an additional condition, color neutrality is guar¬ 
anteed automatically in QCD by the integration over the Ag zero-momentum 
modes. 
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The tadpole diagrams we computed above correspond to the first derivative 
of the effective action at vanishing mean field. 


T° 


5T 


Mg 


A =0 


(5.61) 


In [191] the expectation value of ,4g is computed self-consistently from the Yang- 
Mills equation. 


sr 

Mg 


= 0 . 

A=A 


(5.62) 


This approach is completely equivalent to our approach, which can be seen from 
the Taylor expansion 


5T 

6T 

5 2 r 

Mg 

so 

t-O 

1 

o 

II 

-oo 

l -o 

so 

t-O 

II 


(5.63) 


On the right hand side the first term vanishes, while the second term is essen¬ 
tially the inverse gluon propagator at zero momentum times the expectation 
value of A. 

The tadpole diagram in Fig. 5.2 which we evaluated above corresponds to 
the fermionic part of the color charge density. 


N f 

P(q) =J2^f Ta ^f- 


f= 1 


(5.64) 


One might ask whether the gluons could also contribute to the charge density. 
In Ref. [142] it was argued that the gluonic contribution should vanish, since the 
gluonic part of the charge density 

P a {9 i) = f abcA i F io (5-65) 

contains the chromoelectric field strength, which vanishes in any (super-con¬ 
ducting system. The aim of the next subsection is to corroborate this argument 
by analyzing the tadpole diagram with a gluon loop (Fig. 5.3) [192], 


5.6.4 Gluon self energy and gluon loop tadpole (2SC) 

Let us consider the tadpole diagram in Fig. 5.3, where the rhombus indicates 
a resummed gluon propagator. In [139] integral representations for the vari¬ 
ous components of the gluon self energy in the 2SC phase have been derived. 
(The Feynman diagram for the gluon self energy at leading order looks like the 
second one of Fig. 3.1, with the normal fermion propagators now replaced by 
Nambu-Gor'kov propagators.) It turns out the the gluon self energy in the 2SC 
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Figure 5.3: Tadpole diagram with resummed gluon propagator. 


phase has some off-diagonal components with respect to the color indices, there¬ 
fore the tadpole diagram in Fig. 5.3 does not vanish a priori. The non-vanishing 
components of the gluon self energy in the 2SC phase are 


fin — II22 — IT33, II44 — n 55 — Ilgg — 1177, II45 — —IT54 — IIg7 — —Il7g. 

(5.66) 

In Ref. [139] the gluon self energy is diagonalized by a unitary transformation in 
color space. For our purposes it is more convenient not to perform this transfor¬ 
mation, since we want to keep the structure constants totally antisymmetric. For 
the 44 and 45 components of the self energy analytical results can be obtained 
rather easily, as we will demonstrate in the following. We shall always assume 
that the energy and the momentum of the gluon are much smaller than p, be¬ 
cause if there is any contribution to the tadpole diagrams, it will come from soft 
gluons 5 . In the notation of Ref. [139] we have 


FT 00 - — — n 2 
n 44 — 2 ^ 


/ 


d 3 k 
(27r) 3 


^ (l + eie 2 ki • k 2 ) 

ei,e 2 =± 


x (n°(l - n 2 ) + (1 - 714 ) 712 ) 


1 

Po + e? + e 2 + ip 


1 

po - e? - e 2 + irj 


(5.67) 


and 


-in 


00 

45 


= n 00 = 



d 3 k 

(27t) 3 


^ (1 + e 1 e 2 k 1 • k 2 ) 

ei,e2=± 


x ( 774(1 - 712) - (1 - n°)n 2 ) 


1 

po + e? + £2 + it] 


1 

po - e\ - e 2 + ii ? 


(5.68) 


Here we have used the notations of Ref. [139], 


e* = \J{p - eiki ) 2 + \(j>i \ 2 , 


(5.69) 


rii = 


€j + p- eiki 
2 €i 


(5.70) 


5 For po <f> the gluon self energy is the same as in the normal phase. Since it is then pro¬ 
portional to Sab, it gives no contribution to the tadpole diagram. Furthermore the typical gluon 
momenta relevant for the tadpole diagram would be of the order of the Debye or Meissner masses 
(or even smaller). 
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where (j>i = p ei (p, k t ) is the on-shell gap function. The superscript "0" means 
4>i —> 0. First let us evaluate hull™. Without loss of generality we assume po > 0. 
For po < F the only contribution arises from e\ = e 2 = +. With ki = k + P/ 
k 2 = k, k • p = t, k\ = \J k 2 +p 2 + 2pt, £ = k — p, 4> '■= <t >2 w e find 


2 roo 


d£,{p + £) 2 / dtQ{p 0 -(j))[l + 


x&(p 0 -Ve+^) 


i 

L 2 V / ^Td I 


f + £ + pt 
k\ 


d (po ~ F + h - Vt 2 + 0 2 ) 


+ 


2 V / ^T^ 


<5 (po + F - Ah — VC 2 + </> 2 ) 


(5.71) 


Because of the step function in the second line we have £ <C p for po <C p. 
Therefore it is sufficient to approximate <f> with its value at the Fermi surface, 
which we denote with 4>q- Rewriting the delta functions as 6 (t — ...), we find for 

p <C p 


Imltri = -0(p o - </>o): 


"VPcH’o 


16vrp 

(k + 0 (2p - po + Pt + £ + «) 


d£ (f + £)' 


dt 


'-1 
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«(p + £) 

(« - 0 (2p + PO + pt + g - k) 

«(f + 0 


0(p + Po + £ - «) 0(p - Po - £ + «) S{t - ti) 

0(p + Po ~ £ - «) 0(p - Po + £ + «) <5(i - £ 2 ) 

(5.72) 


with k = an h 

2p (-p 0 - £ + «) + 0o “ P 2 + Po ~ 2 Po« 

1_ 2p(p + 0 

T = 2p (to — £ — *0 + do — P 2 + Po — 2p 0 K 
2 p(p + 0 

We make the substitution £ —» — £ in the term with ip, which leads to 


(5.73) 

(5.74) 


nn Q 2 /VPoWo 

Imn“ = — 0(p o - / _df 0(p + Po - £ - «) ©(p 

iGvrp 

x— [(4p 2 — p 2 + (po + 2£) 2 )(—£ + k) — do(2po + — k)] • 

K 

The ^-integration is now straightforward, and we find 


Po + £ + «) 

(5.75) 


ImFl^ = 0(p o - 0o)0(-Po + y/p 2 + ^o) 

X 2^ \/ p o-^o ( 12 h 2 " 3P 2 + Po - <^o) 
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-0(po - \Jp 2 + 4> t 


g 2 P 2 ^l 


1 + 


2vr {Po~p 2 ) 

-6(f>lPo(Po ~ P 2 ) + 3(Po “ P 2 ) 3 + 0o(P 2 + 3 Po) 


(5.76) 


12/.i 2 (pQ — p 2 ) 2 

This result is a generalization of the one-loop gluon self energy at zero tempera¬ 
ture in the normal phase given in [77] (see also Sec. 3.2). As in the normal phase, 
the leading part of the gluon self energy is of the order g 2 fi 2 (for po. p <C //). In 
order to maintain consistency in the following, we shall keep only these leading 
terms, since contributions from lower powers of p might mix with contributions 
from diagrams with a higher number of loops. Denoting this approximation 
with a tilde, we have 


Imri™ = -&(po - <j>o) 


g 2 p 2 

2n 


©(-Po + 


p 2 _|_ 


Vp 'o ~ 
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+&{p 0 - \lp 2 + <t>l) 


Po 


p 2 J 


(5.77) 


In the following we will refer to this analogon of the HDL approximation as the 
"leading order" approximation, hull™ is an odd function of po, therefore the 
real part can be calculated with the following dispersion relation 6 [167] 
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ReIl£(po,p) = -V 

7T 


duj Imfl™ 
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(5.78) 


^ + Po w-p 0 . 

In contrast to the 11, 22, 33 and 88 components of the gluon self energy, the 
principal value integral can be performed analytically in this case. 


Ren^4 = - 
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7r 2 




arctan 


+© ipo ~ 4o) 


j + vSEl, 
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2(p 2 ~ Po) 


log 


%+P 2 ~P 2 0 


(5.79) 


As a consistency check, we may extract the Debye mass in the normal and in the 
superconducting phase, 

2 2 

— lim lim lim RefiX(py, p) = — ^ , (5.80) 

p^Opo^O 4>o^0 7 T z 

2 2 

- !im lim Refl^(po,p) = (5.81) 

_p—>Opo^O ITT- 

6 In principle it is not sufficient to take only the part of Eq. (5.67) where ei = e 2 = + when 
computing the real part of the gluon self energy. However, as argued in [167], for po,p <C g an¬ 
tiparticles (being always far from the "Fermi surface") only give a constant term in the transverse 
gluon self energy, and this constant is the same in the superconducting phase and in the normal 
phase (at leading order). 
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which are the standard results for Nf = 2 [139]. For the other Lorentz compo¬ 
nents of the self energy one finds following the same steps as above 


9 V 


ImH& = -0(p o - ^o)^2 P' 


@(p-v) l 2p 0 v + 4>l log 
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(5.82) 


(5.83) 


where v = \/po ~ ^o- Again the real parts could be obtained analytically via 
dispersion relations, but the resulting expressions are rather unwieldy and we 
shall refrain from writing them down here. 

The off-diagonal components of the gluon self energy can be evaluated in an 
analogous way. In place of Eq. (5.75) we find 


Imfl 00 = -Q(pq-4>q) 


W'Po-^o 


4yr P J-VpF+1 


d£ ®(p + Po - £ ~ «) ®(p -P 0 + Z.+K) 


x-(</>o-(po + 2£)(k-£)), 


(5.84) 


which gives 


hiill 00 = — 0(p o — \J p 2 + fifyg 2 p 


2 PO%m+P ~Po) 


2vr(p§ - p 2 ) 2 


(5.85) 


We observe that in this one-loop result there is a term which is linear in p, but 
there is no term of order g 2 p 2 . Therefore we have at our order of accuracy 


Imfl 00 = 0. 


(5.86) 


In fact this can be seen rather directly from the (ei = = +)-component of Eq. 

(5.68): if we set d 3 k —» 2np 2 d^dt we find at the order g 2 p 2 that the integrand is 
odd with respect to £ —► — £, t — > —t, and therefore the integral vanishes at this 

order. This argument also holds for the real part, RelT 00 = 0. In the same way it 
can be shown that fl 01 = 0 and IE = 0. 
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Figure 5.4: Tadpole diagram with resummed gluon propagator and three-gluon vertex 
correction. 


This implies that the resummed gluon propagator is diagonal in the color 
indices at leading order. Since the three-gluon vertex is antisymmetric in the 
color indices at tree level, we conclude that the tadpole diagram of Fig. (5.3) 
vanishes at the order of our computation. This means that at this order the 
expectation value of the gluon field, Eq. (5.58), is not changed by the tadpole 
diagram in Fig. (5.3). Therefore the effective chemical potential is determined 
completely by the quark part of the color charge density at the order of our 
computation, whereas the gluonic part is negligible at this order. 

We note that there is no contribution from the tadpole diagram with a ghost 
loop. The reason is that there is no direct coupling between quarks and ghosts in 
the QCD Lagrangian, therefore the ghost self energy vanishes at the order g 2 p 2 , 
and the ghost propagator is proportional to the unit matrix in color space at 
leading order. Furthermore we remark that e.g. in covariant or Coulomb gauge 
in particular the gauge dependent part of the gluon propagator is diagonal in 
the color indices at leading order. Therefore the effective chemical potential is 
gauge independent at this order. 

5.6.5 Gluon vertex correction (2SC) 

We would like to check whether the above result is modified if one replaces 
the tree level gluon vertex with the one-loop vertex correction, as shown in Fig. 
5.4. We remark that the tadpole diagram in Fig. 5.4 corresponds to a higher 
order correction to th e fermionic part of the charge density [Eq. (5.64)]. Fet us 
evaluate the one-loop three-gluon vertex of Fig. 5.4 in the limit where one of the 
gluon momenta approaches zero. After taking the trace with respect to Nambu- 
Gor'kov, color and flavor indices we find 



2^G+(K 1 )'y°G+(K 1 )YG + (K 2 ) - Y G~ G~ {KJY G^ {K 2 ) 

2Yg^(k 1 ) 1 °g^(k 1 )yg-(k 2 ) + yz + (k 1 )Yz-(k 1 )yg^(k 2 ) 


+Y^~(Kih°E+(K 1 )YG+(K 2 ) 


(5.87) 
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-pl-iOu 
L 181 


2 g_ 3 
a/3 


, d 4 A^ m 

=- 


*(27t) 


7^ + (A- 1 )7°G + (A 1 ) 7 "G + (A- 2 ) 


-7^-(A 1 )7°G-(A^ 1 )7"G-(A 2 ) + 7^ + (A^ 1 )7°H~(A 1 )7"G-(A- 2 ) 
-7^-(A- 1 )7 0 H + (A- 1 ) 7 "G + (A 2 ) + 7^ + (A- 1 )7°G + (A^ 1 )7"H-(A^ 2 ) 
-7^-(A- 1 )7°G-(A 1 ) 7 ^ + (A 2 ) - 7 ^-(A: 1 )7 0 S+( j A 1 )7-S-( j A 2 ) 

+7^G + (A 1 )7°H~(A 1 )7^ + (A 2 ) 


(5.88) 


-pl-iOu 

L 484 


9 

V3. 


d 4 K 

i(2ir) 4 


Tr D 


'fG+(K 1 )'fG+(K 1 )'fG+(K 2 ) 


—2 7 m Gq"(A' 1 ) 7 °Gq"(A' 1 ) 7 l/ G + (A' 2 ) - 7 ^G-( j A 1 )7°G-( j A 1 )7"G 0 -(^ 2 ) 
+ 2 7 ^o (AT 1 ) 7 °Gq (K\)YG- (A 2 ) + 7 ^S+( A 1 ) 7 °S- (K^Gq (A 2 ) 

v^-(X 1 )7 0 S+( j A 1 ) 7 "G+(^ 2 


■p/jOfc' 
1 888 


"7 
2<? 3 


(5.89) 


ci 4 A' 


3v/3 7 i(27r)‘ 


rTro 


rG + {I<i h °G + (I<i )YG + {K 2 ) 


- 47 ^G 0 + (A- 1 ) 7 °G 0 + (A' 1 ) 7 ^ 0 + (A 2 ) - 'fG-{K{)'fG-{K 1 )'fG-{K 2 ) 
+4 7 m Gq (AT 1 ) 7 0 Gq (^ 4)7^0 (A 2 ) + 7^H + (A 1 )7°S- {K x )YG~ (A 2 ) 
-7^-(A^ 1 )7°H+(A- 1 )7-G + (A 2 ) - 7^+(A- 1 )7 0 G + (A- 1 )7"H"(A^ 2 ) 
+^-(Kih°G-(K 1 )Y^ + (K 2 ) + 7 /i G - (A' 1 )7°H + (A' 1 ) 7 I 'H~(A" 2 ) 

-7^G + (A 1 )7 0 H-(A 1 )7"H + (A 2 ) 


(5.90) 


where K\ = A" + A and A 2 = AT. Furthermore one finds the relations 


h485 — -h584 — r 687 — _ r 7 gg, Tigi — r 282 — r 383 , 

r 48 4 = r 585 = r 686 = r 787 . (5.9i) 

The propagators G ± and 6'/ in Eqs. (5.87)-(5.90) are defined in Eqs. (5.35) and 
(5.36), and is given by [139] 


=+(K) 


~-(K) 


EE 

h=r,l e=± 

EE 

h=r,l e=± 


</>%(*) 

k“o ~ (fi — ek) 2 — 

mv]* 

ky — (n — ek ) 2 — 



V- h A 


—e 

k ’ 



(5.92) 

(5.93) 


(In the following we will assume again |0®| = 7f j = 7 e [139].) In order to 
evaluate Eqs. (5.87)-(5.90) we need the following Dirac traces. 


Tfi = Tr D 7 M Ak 1 1 7 0 7°Ak 1 7 0 7 I/ Ak 3 ,7 l 


ki 


ki 


k2 


T (2) = Tr D \^V- hl A^°V h2 A^YK: 3 ^ 


(5.94) 

(5.95) 
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One finds 


rrifliy 

1 (3) 

= 

J (4) 

= 

/Tl00 

J (l) - 

^eie2 

II 

E-i 

rpiO 

2(1) 

II 

■'-jCh' 

'e-t 

^eie2 


Tr D 
Tr d 


'rv-h.A-^A^yrvh^ 

7^,T 0 7%A- 7 r A2 A-; 3 


(l + eie 3 ki • k 2 ) , 

= S eie2 (ei k\ + e 3 k 2 ^j , 

5 lj (l - eie 3 ki • k 2 ^ + eie 3 [k\k{ + k 2 k {^ 


(5.96) 

(5.97) 

(5.98) 

(5.99) 
(5.100) 


T (2)' T (3) and T (4) contain a common factor |d/ M /,. 2 , which we do not write in the 
following. Then we find 


^n00 _ rynOO _ rynOO _ /yiOO 

J (l) - 1 (2) ~ 1 (3) ~ 2(4)5 

rpOi _ rpOi _ rpOi _ rpOi 

J (l) - ~ i (2) - 1 (3) - - J (4)> 

rpiO _ rpiO _ rpiO _ rpiO 

~ 1 (2) ~ ~ 1 (3) ~ J (4)> 






(5.101) 

(5.102) 

(5.103) 

(5.104) 


After performing the energy integration one finds for instance 

'■ 3 r d 3 k 


T%g(P,0,~P) = 


w 


f d 6 k 

2 VsJ ^} l + ' 


+ 


eie 2 ki • k 2 

ei,e 2 =± 

2_\ 2 4>\ + - 3£ieisgn(g 2 ) 

(fo - ei ~ 4 + ir l) 2 (po + ei + e 2 + ip) 2 J 

( 


+ 


4 


\Po - ei- e° 2 + ip p 0 + e 4 + e° 2 + ip) ef 


(5.105) 


T°°°(P,0,-.P) = - 


d 3 k 


2y/3 J (2vr) 3 


+ 


^ (l + eie 2 ki-k 
ei ,e2=i 

1 


(po - ei - + ip) 2 (p Q + ei + + ip) 2 


-sgn(6)(2^f + 3£f) + 3^iei 




po - ei - e 2 + ip Po + ei + e 2 + w// 


\ sgn(C 2 )^> : 


(5.106) 


where ki = k + p, k 2 = k, and 6/ = eiki — //. The integrals can be evaluated in a 
similar way as in the previous section. It is easy to see that van¬ 

ishes at the order g 3 p 2 . In general (see also [140]) one finds that T^(P, 0, —P) 
is non-vanishing at the order g 3 p 2 only for those combinations of color indices 
where f a b c is non-vanishing. In the previous section we have seen that the gluon 
propagator is diagonal in the color indices at leading order. Therefore the tad¬ 
pole diagram in Fig. 5.4 also vanishes at this order. 
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5.6.6 Mixing with Nambu-Goldstone bosons (2SC) 

The symmetry breaking pattern of the 2SC phase is SU (3) c x U(1)b —> SU (2) x 
U(1)b [39]. Therefore five massless (would-be) Nambu-Goldstone (NG) bosons 
appear in this phase, which correspond to fluctuations of the diquark conden¬ 
sate. As shown in [136, 137, 138], one finds the following effective action for 
gluons, ghosts and NG bosons after integrating out the quarks. 


T = j d 4 x £ (g ) + ^Tr log (S 1 + 

= J d 4 x £( g) + ^Tr log^ 1 + ^Tr (S^W) - ^TV 

+pTr (S'y lt WS'y v n v S'y p n‘>) + ..., (5.107) 

6 


where the gluonic part of the Lagrangian is e.g. in covariant gauge given by 

Ag) = ~\V Wa + (5.108) 

where Fjt u is the usual field strength of the gluons, c and c are the (anti-) ghost 
fields, and I)' 4 ’ is the covariant derivative in the adjoint representation (see e.g. 
[190]). The quark propagator S reads, in the notation of Ref. [138], 


S = 


and is given by 


W(x,y) = -i 


[G+]" 1 T- \ 1 
[Go]" 1 J ’ 


(5.109) 


(5.110) 


Here aA is the (Lie algebra valued) Maurer-Cartan one-form introduced in [136], 


^ = (id^+gA^T a )V, 


(5.111) 


where 


V = exp 


f 7 i 

Va T a + 


\a =4 


(5.112) 


parametrizes the coset space SU( 3) c x U(1)b/SU(2) x U(1)b [137], B = (1 + 
v/3Tg)/3 is a generator orthogonal to the one of U(1)b, and <p a are the NG 
bosons. We may expand N' in powers of the fields. 


= —^d»<p 8 B + 
V3 


gK - d^ipa - gfabcA^c - |/ab 8 ^V 8 - Tjfabcftbd^lfc 


-\fa8c^d^(p c - x/a&8<£&<9 / V 8 + . . . 
6 6 


T„ 


(5.113) 
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where (p a = tp a for a = 4, 5,6, 7 and (p a = 0 otherwise. 

Let us examine the various terms in the effective action (5.107). First consider 
the term linear in ft, which contains the one-loop gluon tadpole 7 (Fig. 5.2). In 
principle this term also gives contributions to higher n-point functions, since f l 
contains arbitrarily high powers of ip. In particular, there is a contribution to 
the bosonic self energy, but this effect is negligible at leading order since the 
one-loop tadpole [Eq. (5.45)] is only linear in fi. 

The term quadratic in contains the bosonic self energy. As in section 5.6.5 
one finds that at leading order the self energy is diagonal in the color indices. 
There are mixed terms between gluons and NG bosons, which could be elimi¬ 
nated by choosing a suitable t'FIooft gauge [138]. This is not necessary for our 
purposes, but we note that choosing this t'FIooft gauge would not alter our con¬ 
clusion, since also in the t'Hooft gauge of Ref. [138] the gluon, NG and ghost 
propagators are diagonal in the color indices at leading order. The gluon and 
NG propagators in the t'FIooft gauge are given explicitly in Eqs. (54) and (55) of 
Ref. [138]. At leading order the propagators (and in particular the gauge depen¬ 
dent parts of the propagators) are diagonal in the color indices also without the 
unitary transformation in color space that is employed in Ref. [138]. 

The term quadratic in also contains new types of three-boson vertices pro¬ 
portional to fabc, which involve at least one NG boson. 

The term cubic in fl is only non-vanishing at the order y 3 /t 2 for those color 
indices where f abc / 0. In this thesis we do not consider vertices with more than 
three legs, since these would appear only in two-loop tadpole diagrams. 

To summarize, we find at leading order as in the previous subsections that 
the bosonic propagators are diagonal in the color indices, and that three-boson 
vertices of the type (P. 0, —P) are non-zero only for those combinations of 
color indices where f abc / 0. We conclude that the tadpole diagram in Fig. (5.4) 
vanishes at the order of our computation, even if we take into account NG 
bosons. 


5.6.7 Tadpoles with bosonic loops in the CFL phase 

In the case m s = 0 the situation for the diagrams in Figs. (5.3), (5.4) is similar to 
the 2SC phase. In the CFL phase the one-loop gluon self energy is proportional 
to the unit matrix in color space [128] (as in the normal phase). Therefore the 
diagram in Fig. 5.3 vanishes identically. For the gluon vertex correction one 
finds after taking the trace with respect to Nambu-Gor'kov, color and flavor 
indices 


Km( p A-p) = - 


d A K 


Trm if abc 


24 J i( 2tt) 4 1 — c L 7 At G+(A 1 )7°G 2 + (AT)7"G 2 + ( K 2 ) 

+7^r(AA)7 0 G 2 (AT)7"G 2 -(A 2 ) + 7^+(AA)7 0 G+(AT) 7 "G 2 + (A' 2 ) 


7 The ip-tadpole vanishes because the external external momentum of the tadpole diagram is 


zero. 
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+^(K 1 )^Gi(K 1 )YG^(K 2 ) + ^G+(K 1 ) 1 °G+(K 1 )YGi(K 2 ) 
+rGi(K 1 ) 1 °Gi(K 1 )YGi(K 2 ) + fyyi*G+(K 1 ) 1 0 G+(K 1 )'fG+{K 2 ) 
+6^Gi(K 1 ) 1 0 Gi(K 1 )YGi(K 2 ) + ^E+(K\)^ (K\)Y G^ (I< 2 ) 
+rZ 2 (K 1 h°Zt(K 1 )YGt(K 2 ) - 
- 1^2 (Kih°z+(K 1 )'fG^ (k 2 ) + ( a : i ) 7 0 h + (k 1 )yg+ (k 2 ) 

+^E+(K 1 )^(K 1 )YG^(K 2 )+ 1 ^(K 1 ) 1 °Ei(K 1 )YGt(K 2 ) 
+^E+(K 1 ) 1 °Ei(K 1 )YGi(K 2 ) + 2^Ei(K 1 ) 1 °E+(K 1 )YG+(K 2 ) 
+2/H+(^) 7 °5 2 -(ifi)/G 2 -(X 2 ) + /S+^^G+fifO/S,-^) 
+ 7 ^-(^ 1 ) 7 0 G 2 -(X 1 ) 7 ^+(^ 2 ) + 7 ^+(^ 1 ) 7 °G+(^ 1 )rs-(^ 2 ) 
+rZ 2 (K 1 )'y°Gi(K l )YZt(K 2 ) + 2 7 ^+(A' 1 ) 7 °G+(^ 1 ) 7 ‘'~ 2 (K 2 ) 
+2 7 ^ 2 -(K 1 ) 7 0 G 2 (X 1 ) 7 ^+(X 2 )+ 7 ^ 2 (^ 1 ) 7 °Gr(X 1 )rs+(^ 2 ) 
+rs+(^ 1 ) 7 0 G+(X 1 ) 7 "S 2 (^ 2 ) - y t H 2 (A' 1 ) 7 °G 2 (il' 1 ) 7 ^+(K 2 ) 
-rs+(^ 1 ) 7 0 G+(X 1 ) 7 ^ 2 (^ 2 ) + 7 ^r(Xi) 7 0 S+(^ 1 )rS 2 -( J A 2 ) 
+ 7 ^+( j a 1 ) 7 °h 2 ( j a 1 ) 7 ^+(^ 2 ) - (A 2 ) 

-rG 2 + (X 1 ) 7 0 H 2 (X 1 ) 7 ^+(^ 2 ) + 7 ^+(X 1 ) 7 0 H r (^ 1 )rs+( J A 2 ) 
+ 7 %-(A\) 7 0 H+(A 1 )/S 2 (A 2 ) + /G+(A 1 ) 7 0 S 2 -(A- 1 )/S+(il 2 ) 
+ 7 m G 2 (AT 1 ) 7 °H+ (K^Y^i(K 2 ) + 2'fG+(K 1 )'y 0 Et(K 1 )'y v E+(K 2 ) 

+2rG^(K 1 ) 1 0 E+(K 1 )Y^(K 2 ) 

+d abc 7 m G+( j A 1 ) 7 °G+( j A 1 ) 7 i/ G+( j A 2 ) 

-rGr(^ 1 ) 7 0 G 2 (A 1 ) 7 "G 2 (A: 2 )+rG+(^ 1 ) 7 0 G+(A' 1 ) 7 ^+(A 2 ) 
-rG 2 (^ 1 ) 7 0 Gr(A 1 ) 7 ^ 2 ( J A 2 )+ 7 ^G+(^ 1 ) 7 0 G 2 + (A' 1 ) 7 "G+(A 2 ) 
- 7 ^G 2 -( J A 1 ) 7 0 G 2 (^ 1 )rGr(X 2 ) + 6 ^G+(K 1 )' r 0 G+(K 1 )'fG+(K 2 ) 
-6^G 2 (K 2 ) + 7 M S 2 (ATi) 7 0 H^ (-K’ 1 ) 7 I/ G] _ (I< 2 ) 

-'f~i(K 1 )'y 0 ~t(Ki)'fGt(K 2 ) - ^E+iK^E-iK^YG^K,) 
+rZ 2 (K 1 h 0 Et(K 1 )YG+(K 2 ) + rEi(K 1 ) 1 0 E+(K l )YG+(K 2 ) 
- 7 ^+(A 1 ) 7 0 S 2 (A 1 )rG 2 (A 2 ) + 7 ^ 2 (A 1 ) 7 °S+(^ 1 )rG 2 +(^ 2 ) 
-'fE+(K 1 ) 1 °El(K 1 )'fG2(K 2 ) + 2 7 ^ 2 -(A 1 ) 7 °H+(A' 1 ) 7 "G+(A 2 ) 
-2y t H+(^ 1 ) 7 0 H 2 (A 1 ) 7 I 'G 2 (^ 2 )+ 7 ^+(A' 1 ) 7 °G+(^ 1 ) 7 ^ 2 (A 2 ) 
- 7 ^r(A 1 ) 7 0 G 2 (^ 1 ) 7 ^ 2 + (A' 2 ) + y t H+(A 1 ) 7 °G+(A 1 ) 7 ^r(A 2 ) 
-rS 2 (i^ 1 ) 7 °G 2 -( J A 1 ) 7 "S+(^ 2 ) + 2 7 ^ 2 + (A' 1 ) 7 0 G+(^ 1 )rH 2 (A' 2 ) 
-2 7 ^ 2 -( J A 1 ) 7 0 G 2 ( J A 1 ) 7 ^+( J A 2 )+ 7 ^ 2 (^ 1 ) 7 0 Gr( J A 1 )rs+(^ 2 ) 
-^(A'O/G+^b^(A' 2 ) - 7 ^ 2 -(^ 1 ) 7 °G 2 -(^ 1 )rs+(i^ 2 ) 
+ 7 ^+(A 1 ) 7 0 G+(^ 1 ) 7 ^ 2 (A 2 ) + rGr(^ 1 ) 7 0 H+(A 1 ) 7 ^ 2 (A 2 ) 
-'fGt(K 1 )'y°E^(K 1 )'r v E^(K 2 ) - 'fGt(K 1 )'y 0 E+(K 1 )'fE^(K 2 ) 
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+y i Gl(K 1 ) 1 °Zi(K 1 )YZi(K 2 ) + 7^+(A- 1 ) 7 °Hr(K 1 ) 7 ^+(A' 2 ) 
-7^ 2 -( J A 1 )7°H+( J A 1 )7^ 2 -(^ 2 ) + 7^+(X 1 )7°S 2 -(A: 1 )7 i/ S+(^ 2 ) 
-7 /i G 2 (A^)7 0 H+(A^ 1 )7^r(A' 2 ) + 2 7 ^+( J A 1 )7°S 2 ( J A 1 )7^ 2 + (^ 2 ) 


-27^ 2 -(A 1 )7 0 H+(A^ 1 )7^ 2 - 



(5.114) 


where K\ = K + P, A' 2 = K, and [128] 


Gt(K) 

^n(K) 
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ko T (h - ek) 


e=± 


^E 


kl - (fjt - ek) 2 - r {n) (K) 

*U K ) 


kl ~ (fi ~ ek) 2 - 0( n) (K) 


2 A ±e (k)7o, 

A Te (k) 75 . 


(5.115) 

(5.116) 


(As in the 2SC case we have assumed that right-handed and left-handed gap 
functions are real and equal up to a minus sign [128].) Here 4>(i) is the singlet gap 
and <A 2 )(= is the octet gap [128, 134], It is straightforward to show that 

in the three-gluon vertex correction V^(P, 0, —P) only the term proportional to 
fabc is non-vanishing at the order fj 'ir. 

The structure of effective action for gluons and NG bosons is similar to the 
2SC phase [127,131,133,134], Thus we find that all the three-boson vertices are 
proportional to f a i, r at leading order, and all the boson propagators are propor¬ 
tional to 5 a b■ Therefore the tadpole diagram in Fig. 5.4 vanishes also in the CFL 
phase. 


5.7 Some remarks on the specific heat (outlook) 


We would like to conclude this chapter with some remarks on the specific heat 
of color superconducting phases. First let us consider the contribution of a 
fermionic mode to the specific heat. We assume that the imaginary part of the 
fermion self energy is negligible. In a similar way as in Eq. (4.14) one finds then 
the following contribution to the entropy density. 


5f e 


-2 


cfiq 

(27t) 3 



dq 0 dn f (q 0 ) 
2t r dT 


Im log S (c/o + A?w) \ 


(5.117) 


where S is the propagator of the corresponding fermionic mode 8 . Since we ne¬ 
glect the imaginary part of the fermion self energy, the imaginary part of the 
logarithm in Eq. (5.117) gives essentially a step function 0(c/o — u(q)), where 
u(q) is the solution of the dispersion relation .S' 1 (>'((/)• <7) = 0. For the specific 

8 Compared to the notation of chapter 4, we have shifted the energy variable by p, such that 
now the propagator instead of the distribution function depends on p. This notation is much 
more convenient in the context of superconductivity. 
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heat we find then [see Eq. (4.12)] 


C 


•u,ferm. 



(5.118) 


with e(q) = |tc((/)|. For a gapped fermionic mode we have e(q) =y/(q — p) 2 + </> 2 . 
Then one finds that the corresponding contribution to the specific heat is expo¬ 
nentially suppressed, oc exp (—0/T). For an ungapped mode one finds the usual 
Fermi liquid result, C v oc // 2 T, provided that uj(q) is analytic in q — //. 

For a bosonic mode one finds in a similar way (again neglecting the imagi¬ 
nary part of the bosonic self energy) 


r 

'-"l) 


d 3 q dn b (ti(q)) 
^ BT - 


(2vr 


(5.119) 


where Co(q) is the solution of the respective bosonic dispersion law. For massless 
bosons one finds C v oc T 3 , whereas for massive bosons (m S> T) the contribution 
is exponentially suppressed, oc exp(—m/T). 

In the CFF phase all fermionic quasiparticles are gapped. Therefore the 
specific heat is is at low temperature dominated by light bosonic excitations 
[52], leading to a specific heat proportional to T 3 . Other color superconducting 
phases discussed in the literature have also ungapped fermionic quasiparticles, 
leading to a specific heat proportional to ji 2 T [57]. E.g., in the 2SC phase the 
quarks of color 3 are ungapped. Since they have only couplings with massive 
gluons [139], one expects that there will be no anomalous /r T log T contribution 
to the specific heat [57], Also for the TOFF phase the specific heat is linear in T , 
at least in the NJF model approach [193]. Only for the gapless CFF phase the 
specific heat is of the order /i 2 \R7T as a consequence of an (almost) quadratic 
dispersion relation of one of the fermionic quasiparticles [151]. 

These results for the specific heat will receive corrections from the energy- 
momentum dependence and from the non-vanishing imaginary parts of the self 
energies and the gap. These issues certainly deserve further studies, since the 
specific heat is of central importance for the cooling behavior of compact stars. 
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Chapter 6 

Conclusions and outlook 


In this thesis we have investigated properties of cold dense quark matter. The 
initial motivation was the fact that in the core of (some) neutron stars the den¬ 
sity may be high enough for deconfinement of quarks at comparatively small 
temperatures. 

In chapters 3 and 4 we have computed the quark self energy and the specific 
heat in ultradegenerate QCD. At high density and zero temperature the chromo- 
magnetic screening mass vanishes, and quasistatic chromomagnetic fields are 
only dynamically screened. These long-range interactions lead to a non-Fermi- 
liquid behavior at high density and small temperature. In the quark self energy 
the leading result at zero temperature is proportional to g 2 s (E—g) log (geSl^/(E— 
//)) [see Eq. (3.77)]. Correspondingly, the leading term in the interaction part of 
the specific heat is of the order gl s g?T\og{g eS g/T) [see Eq. (4.163)]. We have 
corrected an error in a recent paper [92], in which it was claimed that the lead¬ 
ing term in the interaction part of the specific heat should instead be of the order 
g 2 T 3 log (g e ffn/T). Furthermore we have performed a perturbative expansion of 
the quark self energy and the specific heat in powers of T/ (g e s/j). These expan¬ 
sions contain fractional powers which come from the dynamical screening scale 
q ~ (.(/eff/^Vo)' '• An important application of the results for the specific heat 
and the quark self energy is neutrino emission from ungapped quark matter in 
neutron stars [57], 

Chapter 5 of the thesis has been devoted to color superconductivity. Using 
a general gauge dependence identity we have given a formal proof of gauge 
independence for the fermionic quasiparticle dispersion relations in a color su¬ 
perconductor. As long as the gauge dependence of the quark self energy can be 
neglected, this implies gauge independence of the gap function on the quasipar¬ 
ticle mass shell. The application of gauge dependence identities for gauge the¬ 
ories with spontaneous symmetry breaking at finite temperature has also been 
demonstrated in chapter 2 for the Abelian and a non-Abelian Eliggs model. We 
found furthermore that the spontaneous breaking of global color symmetry in¬ 
duces a non-vanishing expectation value for the gluon field in a color super- 
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conductor. This expectation value acts as an effective chemical potential for the 
color charge. As shown explicitly in [191] this mechanism ensures color neu¬ 
trality of the color superconducting system. The expectation value of the gluon 
field can be computed from the gluon tadpole diagram. We have computed the 
leading order tadpole diagrams, both for the 2SC phase and for the CFL phase, 
where we have also included a small, but non-vanishing strange quark mass in 
the latter case. We have also shown that the expectation value of the gluon field 
is at leading order not modified by one-loop tadpole diagrams with resummed 
gluon or Nambu-Goldstone boson propagators. 

There are several natural continuations of the present work, some of which I 
would like to sketch briefly. As mentioned already in Sec. 4.6 one could use the 
results of chapters 3 and 4 to compute the cooling behavior of a neutron star with 
a normal quark matter component beyond the leading logarithmic accuracy of 
Ref. [57], Still a lot of work has to be done in the field of color superconductiv¬ 
ity. We have computed the gluon tadpole diagrams only for the 2SC and CFL 
phases. A straightforward exercise would be the computation of gluon tadpoles 
for other color superconducting phases, and at finite temperature. It would cer¬ 
tainly be more challenging to compute higher order corrections to the tadpole 
diagram. Probably such a computation will only be possible after higher order 
corrections to the gap have been determined consistently. Another important 
subject is the specific heat of color superconducting phases, which we briefly 
discussed at the end of chapter 5. Also the inhomogeneous color superconduct¬ 
ing (LOFF) phases deserve further studies. In particular it would be interesting 
to compute the neutrino emissivity of neutron stars which contain quark matter 
in the LOFF phase [57]. 

It is fair to say that at present our understanding of cold dense quark matter 
is still very incomplete. As discussed in the Introduction, the major obstacle is 
the fact that lattice simulations at large chemical potential are not possible so far. 
In view of this situation the (semi-)perturbative methods which we have used 
in this thesis provide an extremely valuable tool for an (at least) approximate 
description of the properties of cold dense quark matter. 
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Appendix A 

Matsubara sums 


In this appendix we shall briefly review some standard techniques for the eval¬ 
uation of Matsubara sums. 

For the computation of Feynman diagrams at finite temperature in the imag¬ 
inary time formalism [76, 115] we need the following Matsubara sum in the 
bosonic case, 

J b :=Tj2f(ko = ^inT). (A.l) 

nEZ 

One can rewrite this sum as a contour integral using the residue theorem. As¬ 
suming that /(ho) is regular along the imaginary axis, one gets 

J b = j dko f(k 0 coth (A.l) 

where we have used the relation 

Res V- coth ^1=1. (A.3) 

ko=2mnT 1 1 

The contour C consists of the small circles shown in Fig. A.l, which enclose the 
singularities of the coth at ko = IninT. The contour can be deformed to two 
vertical lines enclosing the imaginary axis, see Fig. A.l. In this way we arrive at 

-i rioo -i rioo-\-e 

Jb = o -7 / dk 0 f(k 0 ) + — dk 0 (f(k 0 ) + f(-k 0 ))n b (k 0 ). (A.4) 

Flere the first term on the right hand side is the vacuum contribution, and the 
second term is the thermal contribution involving the Bose-Einstein distribution, 

Ub ( k 0) = e k 0 /T_ 1 - ( A ' 5 ) 

The sum (A.l) can be evaluated in a slightly different way by considering 
instead of (A.2) the object 

Jb(r) := dk 0 f(k 0 )e koT n b (k 0 ), (A.6) 
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Re&o 


Figure A.l: Integration contours for bosonic Matsubara sums. 


which fulfills lim r _>o <4(r) = ■//,. Since we have for 0 < r < (3 

lim e k ° T n b (ko) = 0, (A.7) 

ko —»d=oo 

we can add arcs at infinity (loosely speaking) which give no contribution, as¬ 
suming that f(k(i) decays sufficiently fast at infinity In this way we obtain a 
new integration contour that pinches the real axis, as shown by the dashed lines 
in Fig. A.l. Taking the limit r —■> 0+ and assuming that 


f(fa ~ irf) = f* (k 0 + irj) (A.8) 

(Schwarz reflection principle), we find 1 

1 f°° 

J b = - I dk 0 n b (k 0 )Imf(k 0 + i7j). (A.9) 

^ J ~o o 

If the function f(ko) fulfills f(-ko) = f(ko), we can rewrite Eq. (A.9) as 

1 [°° 

J b = - dk 0 (1 + 2 n b (k 0 ))Imf(k 0 + irj). (A.10) 

7T Jo 

We refer to the part of Eq. (A.10) which contains n b as the "n^-part", and to the 
other part as the "non-n^-part". 

For fermions one finds in place of (A.4) 


J f ■.= TY J f{ko = 2m{n + \)T) 

nEZ 

i rioo i rioo+e 

= — / dko f (ko) — ——t / dk 0 (f(ko) + f(-ko))n f (ko), (A.ll) 

_ J _ j - x _Z7T2- J _2oo+e 

1 The point ko = 0 should be excluded from the integration path. 
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and in place of (A.9) 


1 

J f = -/ dk 0 nf(ko)lmf(k 0 + irj), 

^ J — OO 

where rif(ko) is the Fermi-Dirac distribution, 

n S (*») = e t./T + l - 

At finite chemical potential one has to replace nf(ko) with rif ( k (i — 


(A.12) 


(A.13) 
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Appendix B 

HDL/HTL gluon spectral 
densities 


In this appendix we review the spectral densities of the gluon propagator in the 
HDL/ HTL approximation, as given in [76]. From Eq. (3.37) one finds for the 
transverse spectral density 

^Pt{qo, q) = Z T (q) [<%o - v T (q)) ~ %o + u T (q))] + Pr(qo, q). (B.l) 

Here the dispersion law i or{q) is determined from the pole of the propagator, 

co T (q) 2 ~q 2 ~ Ren™ L/HTL (ca T (g), q) = 0, (B.2) 

and the residue Zr(q) is given by 

Wt(4 - q 2 ) 


Z T (q) = 


— (ojj, — q 2 ) 2 


(B.3) 


with the plasma frequency 


LUi 


= V / V3' 


m, 


(B.4) 


where m is defined in Eq. (3.32). For the cut contribution Briqo- q) one finds 

Pr(qo, q) = m 2 x( 1 - x 2 )9( 1 - x 2 )/2 

2 X(1~X 2 ) 


q 2 (x 2 — 1) — m 2 + 

2 4 2 (! -x 2 ) 21 

+7T m x --- 


log 


X + 1 


x — 1 


(B.5) 


where x = qo/q. For the longitudinal spectral density one finds from Eq. (3.38) 

7^PL(qo, q) = Z L (q) [% 0 - u L (q)) ~ S(q 0 + u) L (q))] + (3 L {qo, q). (B.6) 
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Figure B.l: The dispersion laws wx(q) (upper curve) and wi,(q) (lower curve) for lo p = 
1. (Dotted line: light cone). 


Zt Zl 



Figure B.l: The residues Zx(q ) and Zi(q)for w p = 1. (Dotted line in the left panel: -Z, 
dotted line in the right panel: ffj). 


Again the dispersion law u>L(q ) is determined from the pole of the propagator 
[see Eqs. (3.9), (3.38)], 


q 2 + Ren™ L/HTL (a; L (g),g) = 0, 

and the residue Zjfq) is given by 

u l (w 2 l - q 2 ) 

L[q> q 2 {q 2 + 

For the cut contribution q) one finds 


PL(qo,q) = m 

2 x0( 1 — x 2 ) 


X 

(q 2 + 2 m 2 f 1 - | log 

x + 1 

x — 1 


V V 2 


(B.7) 


(B.8) 


n -1 


+ 7r 2 m 4 x 2 


(B.9) 


with x = go/ q- 

Fig. B.l shows the dispersion laws ur{q) and wrf q), as determined from Eqs. 
(B.2) and (B.7). Fig. B.2 shows the residues Zx(q) and Zjf q) which are given by 
Eqs. (B.3) and (B.8). 
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Appendix C 

Matrix identities 


In order to be complete, we give here the proof for two matrix identities which 
are needed in the main text. 

Proposition 1: For the variation of the determinant of an invertible matrix 
M. one has 


5 det M = (det M) Tr [<f log M] = (det M) Tr[M _1 (5M]. (C.l) 


Proof: Using 


we find 

S det M 


det M = exp (TV log M) 

(det M) 5 (TV log M) = (det M) Tr [8 log M]. 


log M can be expanded in a series, for instance about the unit matrix. Under the 
trace cyclic permutations of a product of matrices are allowed. Therefore one 
has 


TV [5 log M] = TrpVrtM], 


□ 


Proposition 2: For any n x n matrices A, B, C, D (with D invertible) one has 


det 


A B 
C D 


= det(DA - DBD~ i C) 


(C.2) 


Proof: With the definition 


E := (A — BD _1 C) -1 


we can write the inverse matrix as 

/ A BV 1 / E EBD 1 \ 

V C D ) “V D _1 CE B _1 AEBD 1 ) ' 


A short calculation using this relation shows that 

-l 


det 


DE’ 1 

D~ 1 C 


A 

C 


B 

D 


= det 


D 

0 


-DBD 1 
D 1 


= 1 . 
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where 0 and 1 denote the n X n zero and unit matrices, respectively. Therefore 
we have 


det 


A B 
C D 


det ( ° E 1( ^ J ) = det (DE- 1 ) . 


□ 
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